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Abstract 

The Green-function technique, termed the irreducible Green func- 
tions (IGF) method, that is a certain reformulation of the equation-of 
motion method for double-time temperature dependent Green func- 
tions (GFs) is presented. This method was developed to overcome 
some ambiguities in terminating the hierarchy of the equations of mo- 
tion of double-time Green functions and to give a workable technique 
to systematic way of decoupling. The approach provides a practical 
method for description of the many-body quasi-particle dynamics of 
correlated systems on a lattice with complex spectra. Moreover, it 
provides a very compact and self-consistent way of taking into account 
the damping effects and finite lifetimes of quasi-particles due to inelas- 
tic collisions. In addition, it correctly defines the Generalized Mean 
Fields (GMF), that determine elastic scattering renormalizations and 
, in general, are not functionals of the mean particle densities only. The 
purpose of this article is to present the foundations of the IGF method. 
The technical details and examples are given as well. Although some 
space is devoted to the formal structure of the method, the emphasis 
is on its utility. Applications to the lattice fermion models such as 
Hubbard/ Anderson models and to the Heisenberg ferro- and antifer- 
romagnet, which manifest the operational ability of the method are 
given. It is shown that the IGF method provides a powerful tool for 
the construction of essentially new dynamical solutions for strongly in- 
teracting many-particle systems with complex spectra. 
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1 Introduction 



The basic problems of field theory and statistical mechanics are much similar 
in many aspects, especially, when we use the method of second quantiza- 
tion and Green functions [||. In both the cases, we are dealing with systems 
possessing a large number of degrees of freedom (the energy spectrum is 
practically a continuous one) and with averages of quantum mechanical op- 
erators IQ]. In quantum field theory, we mostly consider averages over the 
ground state, while in statistical mechanics, we consider finite temperatures 
(ensemble averages) as well as ground-state averages. Great advances have 
been made during the last decades in statistical physics and condensed mat- 
ter theory through the use of methods of quantum field theory @ - [^]- 
It was widely recognized that a successful approximation for determining 
exited states is based on the quasi-particle concept and the Green func- 
tion method. For example, the study of highly correlated electron systems 
has attracted much attention recently @ - @, especially after discovery of 
copper oxide superconductors, a new class of heavy fermions 0, and low- 
dimensional compounds |^] , |^ . Although much work for strongly correlated 
systems has been performed during the last years, it is worthy to remind 
that the investigation of excitations in many-body systems has been one of 
the most important and interesting subjects for the last few decades. 
The quantum field theoretical techniques have been widely applied to sta- 
tistical treatment of a large number of interacting particles. Many-body 
calculations are often done for model many-particle systems by using a per- 
turbation expansion. The basic procedure in many-body theory [^] is to 
find a suitable unperturbed Hamiltonian and then to take into account a 
small perturbation operator. This procedure that works well for weakly 
interacting systems needs a special reformulation for many-body systems 
with complex spectra and strong interaction. For many practically interest- 
ing cases {e.g. in quantum chemistry problems ), the standard schemes of 
perturbation expansion must be reformulated greatly |11| - |15|. Moreover, 
many-body systems on a lattice have their own specific features and in some 
important aspects differ greatly from continuous systems. 
In this review that is largely pedagogical we are primarily dealing with the 
spectra of elementary excitations to learn about quasi-particle many-body 
dynamics of interacting systems on a lattice. Our analysis is based on the 
equation-of-motion approach, the derivation of the exact representation of 
the Dyson equation and construction of an approximate scheme of calcula- 
tions in a self-consistent way. In this review only some topics in the field 
are discussed. The emphasis is on the methods rather than on a detailed 
comparison with the experimental results. We attempt to prove that the ap- 
proach we suggest produces a more advanced physical picture of the problem 
of the quasi-particle many-body dynamics. 

The most characteristic feature of the recent advancement in the basic re- 



3 



search on electronic properties of solids is the development of variety of the 
new classes of materials with unusual properties: high- Tc superconductors, 
heavy fermion compounds, complex oxides, diluted magnetic semiconduc- 
tors, perovskite manganites, etc. Contrary to simple metals, where the fun- 
damentals are very well known and the electrons can be represented so that 
they weakly interact with each other, in these materials, the electrons in- 
teract strongly, and moreover their spectra are complicated, i.e. have many 
branches. This gives rise to interesting phenomena such as magnetism, 
metal-insulator transition in oxides, heavy fermions, colossal negative mag- 
netoresistance in manganites, etc., but the understanding of what is going 
on is in many cases only partial. 

The subject of the present paper is a microscopic many-body theory of 
strongly correlated electron models. A principle importance of these stud- 
ies is concerned with a fundamental problem of electronic solid state theory, 
namely with the tendency of 3d(4d) electrons in transition metal compounds 
and 4f(5f) electrons in rare-earth metal compounds and alloys to exhibit 
both the localized and delocalized (itinerant) behaviour. Interesting elec- 
tronic and magnetic properties of these substances are intimately related to 
this dual behaviour of electrons 1 17 ]-[ 19 1. 

The problem of adequate description of strongly correlated electron systems 
has been studied intensively during the last decade |2C],pl|], especially in 
context of the physics of magnetism, heavy fermions and high- Tc supercon- 
ductivity ||7| . The understanding of the true nature of electronic states and 
their quasi-particle dynamics is one of the central topics of the current ex- 
perimental and theoretical studies in the field. A plenty of experimental 
and theoretical results show that this many-body quasi-particle dynamics 
is quite nontrivial. A vast amount of theoretical searches for a suitable de- 
scription of strongly correlated fermion systems deal with simplified model 
Hamiltonians. These include, as workable patterns, the single-impurity An- 
derson model (SIAM) and Hubbard model. In spite of certain drawbacks, 
these models exhibit the key physical feature: the competition and interplay 
between kinetic energy (itinerant) and potential energy (localized) effects. 
A fully consistent theory of quasi-particle dynamics of both the models is 
believed to be crucially important for a deeper understanding of the true 
nature of electronic states in the above-mentioned class of materials. In 
spite of experimental and theoretical achievements, it remains still much to 



be understood concerning such systems |T^,|22]. 

Recent theoretical investigations of strongly correlated systems have brought 
forth a significant variety of the approaches to solve these controversial prob- 
lems. There is an important aspect of the problem under consideration, 
namely, how to take adequately into account the lattice (quasi-localized) 
character of charge carriers, contrary to simplified theories of the type of a 
weakly interacting electron gas. To match such a trend, we need to develop 
a systematic theory of correlated systems, to describe, from the first princi- 
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pies of the condensed matter theory and statistical mechanics, the physical 
properties of this class of materials. 

In previous papers, we set up the practical technique of the method of the 
irreducible Green functions (IGF) -[!!]• This IGF method allows one 
to describe quasi-particle spectra with damping for systems with complex 
spectra and strong correlation in a very general and natural way. This 
scheme differs from the traditional methods of decoupling or terminating 
an infinite chain of the equations and permits one to construct the relevant 
dynamic solutions in a self-consistent way on the level of the Dyson equation 
without decoupling the chain of the equations of motion for the double-time 
temperature Green functions. The essence of our consideration of dynamic 
properties of many-body system with strong interaction is related closely 
with the field theoretical approach, and we use the advantage of the Green- 
function language and the Dyson equation. It is possible to say that our 
method emphasizes the fundamental and central role of the Dyson equation 
for the single-particle dynamics of many-body systems at finite tempera- 
ture. This approach has been suggested as essential for various many-body 
systems, and we believe that it bears the real physics of interacting many- 
particle interacting systems [^], |25]. 

It is the purpose of the present paper to introduce the concepts of irre- 
ducible Green functions (or irreducible operators) and Generalized Mean 
Fields (GMF ) in a simple and coherent fashion to assess the validity of 
quasi-particle description and mean field theory. The irreducible Green func- 
tion method is a reformulation of the equation-of-motion approach for the 
double-time thermal GFs, aimed of operating with the correct functional 
structure of the required solutions. In this sense, it has all advantages and 
shortcomings of the Green-function method in comparison, say, with the 
functional integration technique, that, in turn, has also its own advantages 
and shortcomings. The usefulness of one or another method depends on the 
problem we are trying to solve. For the calculation of quasi-particle spec- 
tra, the Green- function method is the best. The irreducible-Green- function 
method adds to this statement: "for the calculation of the quasi-particle 
spectra with damping" and gives a workable recipe how to do this in a self- 
consistent way. 

The distinction between elastic and inelastic scattering effects is a funda- 
mental one in the physics of many-body systems, and it is also reflected 
in a number of other ways than in the mean-field and finite lifetimes. The 
present review attempts to offer a balanced view of quasi-particle interaction 
effects in terms of division into elastic- and inelastic-scattering characteris- 
tics. For this aim, in the present paper, we discuss the background of the 
IGF approach more thoroughly. To demonstrate the general analysis, we 
consider here the calculation of quasi-particle spectra and their damping 
within various types of correlated electron models to extend the applicabil- 
ity of the general formalism and show flexibility and practical usage of the 
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IGF method. 



2 Varieties of Green Functions 

It is appropriate to remind the ideas underlying the Green- function method, 
and to discuss briefly why they are particularly useful in the study of inter- 
acting many-particle systems. 

The Green functions of potential theory were introduced to flnd the fleld 
which is produced by a source distribution {e.g. the electromagnetic field 
which is produced by current and charge distribution). The Green func- 
tions in field theory are the so-called propagators which describe the tempo- 
ral development of quantized fields, in its particle aspect, as was shown by 



Schwinger in his seminal works |35] - |37]. The idea of the Green function 
method is contained in the observation that it is not necessary to attempt to 
calculate all the wave functions and energy levels of a system. Instead, it is 
more instructive to study the way in which it responds to simple perturba- 
tions, for example, by adding or removing particles, or by applying external 
fields. 

There is a variety of Green functions |^ and there are Green functions for 
one particle, two particles..., n particles. A considerable progress in study- 
ing the spectra of elementary excitations and thermodynamic properties of 
many-body systems has been for most part due to the development of the 
temperature-dependent Green-functions methods. 



2.1 Temperature Green Functions 

The temperature dependent Green functions were introduced by Matsub- 
ara [^]. He considered a many-particle system with the Hamiltonian 

(1) H = Ho + V 

and observed a remarkable similarity that exists between the evaluation of 
the grand partition function of the system and the vacuum expectation of 
the so-called S-matrix in quantum field theory 

(2) Z = Tr exp[(/iAf - Ho)(3]S{py, S{(5) = 1- T y(r)5(r)dr 

Jo 

where /3 = {kT)~^. In essence, Matsubara observed and exploited, to great 
advantage, formal similarities between the statistical operator exp(—fiH) 
and the quantum-mechanical time-evolution operator exp(iHt). As a result, 
he introduced thermal ( temperature-dependent ) Green functions which we 
call now the Matsubara Green functions. 

We note that the thermodynamic perturbation theory has been invented by 
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Peierls j39U . For the free energy of a weakly interacting system he derived 
the expansion up to second order in perturbation: 

(3) F = Fo + + -^Y: ^nnPn + f (E ^nnPn) " 

n m.n n m n n 

where p„ = exp[/3(Fo - E^)] and exp(-/3Fo) = X^n exp(-/3F°). By using 
the expansion of S{(3) up to second order 

(4) S{p) = 1- t V(T)dT + f dn r dT2V{T{)V(T2) + ... 

Jo Jo Jo 

and rearranging the terms in the expression for Z, it can be shown that the 
Peierls result for the thermodynamic potential Q can be reproduced by the 
Matsubara technique (for a canonical ensemble). 

Though the use of Green functions is related traditionally with the perturba- 



tion theory through the use of diagram techniques, in paper |35| a prophetic 
remark has been made: 

"... it is desirable to avoid founding the formal theory of the 
Green functions on the restricted basis provided by the assump- 
tion of expandability in powers of coupling constants" . 

Since the most important aspect of the many-body theory is the necessity of 
taking properly into account the interaction between particles, that changes 
( sometimes drastically ) the behaviour of non-interacting particles, this re- 
mark of Schwinger is still extremely actual and important. 
Since that time, a great deal of work has been done, and many different vari- 
ants of the Green functions have been proposed for studies of equilibrium 
and non-equilibrium properties of many-particle systems. We can mention, 
in particular, the methods of Martin and Schwinger |3^ and of Kadanoff and 
Baym [^]. Martin and Schwinger formulated the GF theory not in terms of 
conventional diagrammatic techniques, but in terms of functional-derivative 
techniques that reduces the many-body problem directly to the solution of 
a coupled set of nonlinear integral equations (see also[H]). The approach 
of Kadanoff and Baym establishes general rules for obtaining approxima- 
tions which preserve the conservation laws ( sometimes called conserving 
approximations Q ). As many transport coefficients are related to conser- 
vation laws, one should take care of it when calculating the two-particle and 
one-particle Green functions [^]. The random-phase approximation, that is 
an essential point of the whole Kadanoff-Baym method, does this and so 
preserves the appropriate conservation laws. It should be noted, however, 
that the Martin-Schwinger and Kadanoff-Baym methods in their initial form 
were formulated for treating the continuum models and should be adapted 
to study lattice models, as well. 



However, as was claimed by Matsubara in his subsequent paper ||4^, the 
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most convenient way to describe the equilibrium average of any observable 
or time-dependent response of a system to external disturbances is to ex- 
press them in terms of a set of the double-time, or Bogoliubov-Tyablikov, 
Green functions. 

The aim of the present paper is to suggest and justify that an approach , 



the irreducible Green functions (IGF) method [43|, |24|, that is in essence 
a suitable reformulation of an equation-of-motion approach for the double- 
time temperature-dependent Green functions provides an effective and self- 
consistent scheme for description the many-body quasi-particle dynamics of 
strongly interacting many-particle systems with complex spectra. This IGF 
method provides some systematization of approximations and removes (at 
least partially) the difficulties usually encountered in the termination of the 
hierarchy of equations of motion for the GF. 

2.2 Double-time Green Functions 

In this Section, we briefly review the double-time temperature-dependent 
Green functions . 

The double-time temperature-dependent Green functions were introduced 
by Bogoliubov and Tyablikov |44| and reviewed by Zubarev (4|] and Tyab- 



likov |6|. 



Consider a many-particle system with the time-independent Hamiltonian 
a = H — fiN; ^ is the chemical potential, is the operator of the total 
number of particles, and we have chosen our units so that fi = \. Let A{t) 
and B{t') be some operators . The time development of these operators in 
the Heisenberg representation is given by: 

(5) A{t) = exp(iHt)A(0) exp(-iHt) 

We define three types of Green functions, the retarded, advanced, and causal 
Green functions: 

(6) C =« A{t),B{t') >y= -i9{t - t') < [A{t)B{t% >,!] = ±1. 

(7) =« A{t),B{t') »"= ie{t' -t)< [A{t)B{t% >,ri = ±1. 

(8) =« Ait),B{t') »"= iT < A{t)B{t') >= 
ie{t - t') < A{t)B{t') > +7]ie{t' -t) < B{t')A{t) >,r/ = ±1. 

where < ... > is the average over a grand canonical ensemble, 9{t) is a step 
function, and square brackets represent the commutator or anticommutator 

(9) [A,B]r^ = AB -r]BA 

Differentiating a Green function with respect to one of the arguments, for ex- 
ample, the first argument, we can obtain the equation (equation-of-motion) 
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describing the development of this function with time 

id/dtC^it^t') = 6{t-t') < [A,B]n > + « [A, H]{t), B{t') »";a = r,a,c 
(10) 

Since this differential equation contains an inhomogeneous term with 6-type 
factors, we are dealing formally with the equation similar to the usual one 
for the Green function [Q| and for this reason, we use the term the Green 
function. We note that the equation of motion is of the same functional 
form for all the three types of Green functions ( i.e. retarded, advanced, 
and causal ). However, the boundary conditions for t are different for the 
retarded, advanced, and causal functions |44]. 

The next differentiation gives an infinite chain of coupled equations of motion 

(11) (i)"(i"/(irG(t, t') = 



Y,{i)'"'''d''-^/dt"-Hit-t') < [[...[A,H]...H],B]^ > 
+ « [[...[A,H]...H]{t),B{t')» 



To solve the differential equation-of-motion, we should consider the Fourier 
time transforms of the Green functions: 

(12) GAB{t - t') = (27r)-i / duoGAsi^) exp[-iu;(t - t% 

J oo 

roo 

(13) Gab{oj)=«A\B»^= dtGAB{t)exp{iu;t), 

J oo 

By inserting into (0) and (|l|), we obtain 

(14) ioGABiio) =< [A,B]r, > + « [A,H]\B »^; 

n 

(15) io"^GABiu^) = E'^""' < [[.■.[A,H]...H],B]r, > 

+ « [[...[A,H]...H]\B»^ 

n 

It is often convenient to differentiate of the Green function with respect to 
the second time t' . In terms of Fourier time transforms, the corresponding 
equations of motion read 

(16) - ivGABiuj) = -<[A, B]r, > + « A\[B,H] »^; 

n 

i-lTu^^GABiu^) = -^(-1)"-'=^"-'= < [A, [...[B,H]...H]]r, > 

(17) + «A\[...[B,H]...H]»^ 
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It is rather difficult problem to solve the infinite chain of coupled equa- 
tions of motion (|l^) and (0) . It is well established now that the usefulness 
of the retarded and advanced Green functions is deeply related with the 



dispersion relations |44], that provide the boundary conditions in the form 



of spectral representations of the Green functions. 
2.3 Spectral Representations 

The GFs are linear combinations of the time correlation functions: 

I r+oo 



\ r+oo 

{IS^ABit - t') =< A{t)B{t') >= — / dujexplMt - t')]AABiuj) 

\ f + OO 

(19) FBA(t' - t) =< B{t')A{t) >=— duj exp[iu;(t' - t)]ABAiuj) 

2TT J-oo 

Here, the Fourier transforms Aab{^) and Aba{^) are of the form 

(20) ABAiio) = 
Q-^27rJ2eM-PEn){4^iBiJm){i^lAi^n)S{En -Em-uj) 

m,n 

(21) Aab = exp{-puj)ABA{-u^) 



The expressions ( [20| ) and ( ^l| ) are spectral representations of the correspond- 
ing time correlation functions. The quantities Aab and Aba are spectral 
densities or spectral weight functions. 
It is convenient to define 

(22) FBAiO) =< B{t)A{t) >= — / dioA{u;) 

1 r+oo 

(23) Fab{0) =< A{t)B{t) >= — dujexp{Puj)A{u;) 

2tT J-oo 



Then, the spectral representations of the Green functions can be expressed 
in the form 

(24) G'{u;)=«A\B»l= 

[exp{Pu;')-v]A{u;') 



2n J-oo to — to' + ie 
(25) Ciu) =« A\B »l= 

1 i' + OO .1 

[exp{f3u:')-r^]A{u:') 



The most important practical consequence of spectral representations for the 
retarded and advanced GFs is the so-called spectral theorem. The spectral 
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theorem can be written as 



(26) < B{t')A{t) >= 
/ duj exp[iuj{t — t')][exp{/3uj) — r]] ImGAB{^ + i^) 

(27) < A{t)B{t') >= 
/ exp(/3a;) exp[iu;(f — t')][exp(/?a;) — ry]~ ImGABi'-^ + i^) 

Expressions (26) and ([27|) are of fundamental importance. They directly re- 



late the statistical averages with the Fourier transforms of the corresponding 
GFs. The problem of evaluating the latter is thus reduced to finding their 
Fourier transforms , providing the practical usefulness of the Green functions 



technique |45|, 



3 Irreducible Green Functions Method 

In this Section, we discuss the main ideas of the IGF approach which allows 
one to describe completely quasi-particle spectra with damping in a very 
natural way. 



We reformulated the two-time GF method |43], [24| to the form, which is 



especially adjusted [pi], |43] for correlated fermion systems on a lattice and 



systems with complex spectra [^6|,[27|. A similar method was proposed 



in paper |47] for Bose systems ( anharmonic phonons and spin dynamics of 
Heisenberg ferromagnet ) . The very important concept of the whole method 
is the Generalized Mean Field (GMF), as it was formulated in ref. |24]. 



These GMFs have a complicated structure for the strongly correlated case 
and complex spectra and are not reduced to the functional of mean densi- 
ties of the electrons or spins when one calculates excitation spectra at finite 
temperatures. 



3.1 Outline of IGF Method 

To clarify the foregoing, let us consider a retarded GF of the form |4(:] 

(28) =« A{t),A\t') »= -i6{t - t') < [A{t)A^{t% >,r? = ±1 

As an introduction to the concept of IGFs, let us describe the main ideas of 
this approach in a symbolic and simplified form. To calculate the retarded 
GF G{t — t'), let us write down the equation of motion for it: 

(29) ujG{lj) =< [A, AX > + « [A, H]^ \ A^ »^ 



The essence of the method is as follows [24]: 

It is based on the notion of the "IRREDUCIBLE" parts of GFs (or the 
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irreducible parts of the operators, A and , out of which the GF is con- 
structed) in terms of which it is possible, without recourse to a truncation 
of the hierarchy of equations for the GFs, to write down the exact Dyson 
equation and to obtain an exact analytic representation for the self-energy 
operator. By definition, we introduce the irreducible part (ir) of the GF 

(30) « [A,H]^\A^ »=« [A,H]_-zA\A^ » 
The unknown constant z is defined by the condition (or constraint) 

(31) <[[A,Ht\AX>=0 

which is an analogue of the orthogonality condition in the Mori formalism ( 
see ref.|^^). From the condition (^) one can find: 

_ <[[A,HU,A^,> _ Ml 
^ ' <[A,A%> -Mo 

Here Mq and Mi are the zeroth and first order moments of the spectral 
density. Therefore, the irreducible GFs are defined so that they cannot be 
reduced to the lower-order ones by any kind of decoupling. It is worth not- 
ing that the term "irreducible" in a group theory means a representation of 
a symmetry operation that cannot be expressed in terms of lower dimen- 
sional representations. Irreducible (or connected ) correlation functions are 
known in statistical mechanics (c/.|^l|). In the diagrammatic approach, the 
irreducible vertices are defined as graphs that do not contain inner parts 
connected by the G^-line. With the aid of the definition (30) these concepts 
are translated into the language of retarded and advanced GFs. This proce- 
dure extracts all relevant (for the problem under consideration) mean-field 
contributions and puts them into the generalized mean-field GF which is 
defined here as 

(33) O0,,).ll^ 

To calculate the IGF « [A, H]-{t), A^ {t') » in (|9|), we have to 

write the equation of motion for it after differentiation with respect to the 
second time variable t'. It should be noted that the trick of two-time differ- 
entiation with respect to the first time t and second time t' (in one equation 
of motion) was introduced for the first time by Tserkovnikov [p9| . 
The condition of orthogonality (31) removes the inhomogeneous term from 
this equation and is a very crucial point of the whole approach. If one in- 
troduces the irreducible part for the right-hand side operator as discussed 
above for the "left" operator, the equation of motion (|29|) can be exactly 
rewritten in the following form 

(34) G = G° + G°PG° 
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The scattering operator P is given by 

(35) P={Mo)-\ « [A,HU[A\H]^ »(-))(Mo)-i 

The structure of equation (34) enables us to determine the self-energy op- 
erator M, by analogy with the diagram technique 

(36) P = M + MG^P 



From the definition ( |36D it follows that the self-energy operator M is defined 
as a proper (in the diagrammatic language, "connected") part of the scat- 
tering operator M = {Py. As a result, we obtain the exact Dyson equation 
for the thermodynamic double-time Green functions: 

(37) G = G° + G°MG 

The difference between P and M can be regarded as two different solutions 
of two integral equations (|3^ and (^). But from the Dyson equation (|37| ) 
only the full GF is seen to be expressed as a formal solution of the form 

(38) G = [(G°)-i - M]~i 

Equation (38) can be regarded as an alternative form of the Dyson equation 
(37) and the definition of M provided that the generalized mean- field GF 



G^ is specified. On the contrary , for the scattering operator P, instead of 
property G^G~^ + G^M = 1, one has the property 

(G^)-^ -G-^ =PG°G-^ 

Thus, the very functional form of the formal solution (38) determines the 
difference between P and M precisely. 

Thus, by introducing irreducible parts of GF (or irreducible parts of the 
operators, out of which the GF is constructed) the equation of motion (|29| ) 
for the GF can exactly be ( but using orthogonality constraint (|3l|)) trans- 
formed into the Dyson equation for the double-time thermal GF (|37|). This 
result is very remarkable , because the traditional form of the GF method 
does not include this point. Notice that all quantities thus considered are 
exact. Approximations can be generated not by truncating the set of cou- 
pled equations of motions but by a specific approximation of the functional 
form of the mass operator M within a self-consistent scheme, expressing M 
in terms of initial GF 

M ^ F[G] 

Different approximations are relevant to different physical situations. 



The projection operator technique |5C] has essentially the same philoso- 



phy, but with using the constraint ( plf ) in our approach we emphasize the 
fundamental and central role of the Dyson equation for the calculation of 
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single-particle properties of many-body systems. The problem of reducing 
the whole hierarchy of equations involving higher-order GFs by a coupled 
nonlinear set of integro-differential equations connecting the single-particle 
GF to the self-energy operator is rather nontrivial ( c/.|^l|]). A characteristic 
feature of these equations is that, besides the single-particle GF, they in- 
volve also higher-order GF. The irreducible counterparts of the GFs, vertex 
functions, etc, serve to identify correctly the self-energy as 

M = Go ^ - G-^ 



The integral form of Dyson equation gives M the physical meaning 
of a nonlocal and energy-dependent effective single-particle potential. This 
meaning can be verified for the exact self-energy through the diagrammatic 
expansion for the causal GF. 

It is important to note that for the retarded and advanced GFs, the notion 



of the proper part M = (P)^ is symbolic in nature |24]. In a certain sense 



it is possible to say that it is defined here by analogy with the irreducible 



many-particle T-matrix[41|. Furthermore, by analogy with the diagram- 
matic technique, we can also introduce the proper part defined as a solution 
to the integral equation p^). These analogues allow us to understand better 
the formal structure of the Dyson equation for the double-time thermal GF 
but only in a symbolic form . However, because of the identical form of the 
equations for GFs for all three types ( advanced, retarded, and causal ), we 
can convert in each stage of calculations to causal GF and, thereby, confirm 



the substantiated nature of definition (36)! We therefore should speak of an 



analogy of the Dyson equation. Hereafter, we drop this stipulating, since 
it does not cause any misunderstanding. In a sense, the IGF method is a 
variant of the Gram-Schmidt orthogonalization procedure (see Appendix A 
)• 

It should be emphasized that the scheme presented above gives just a gen- 
eral idea of the IGF method. A more exact explanation why one should not 
introduce the approximation already in P, instead of having to work out M, 
is given below when working out the application of the method to specific 
problems. 

The general philosophy of the IGF method is in the separation and identi- 
fication of elastic scattering effects and inelastic ones. This latter point is 
quite often underestimated, and both effects are mixed. However, as far as 
the right definition of quasi-particle damping is concerned, the separation of 
elastic and inelastic scattering processes is believed to be crucially important 
for many-body systems with complicated spectra and strong interaction. 
From a technical point of view, the elastic GMF renormalizations can exhibit 
quite a nontrivial structure. To obtain this structure correctly, one should 
construct the full GF from the complete algebra of relevant operators and 
develop a special projection procedure for higher-order GFs in accordance 
with a given algebra. Then the natural question arises how to select the 
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relevant set of operators {Ai,A2, ■■■An} , describing the "relevant degrees of 
freedom". The above consideration suggests an intuitive and heuristic way 
to the suitable procedure as arising from an infinite chain of equations of 



motion (14). Let us consider the column 



/Ai\ 
A2 



\An) 



where 



Ai = A, A2 = [A,Hl A^ = [[A,HlHl^^^An = [l^^[A,H]..^H] 



Then the most general possible Green function can be expressed as a matrix 



(A^\ 
A2 



\An) 



(4 4 



» 



This generalized Green function describes the one-, two- and n-particle 
dynamics. The equation of motion for it includes, as a particular case, the 
Dyson equation for single-particle Green function, the Bethe-Salpeter equa- 
tion, which is the equation of motion for the two-particle Green function 
and which is an analogue of the Dyson equation, eic. The corresponding 
reduced equations should be extracted from the equation of motion for the 
generalized GF with the aid of the special techniques such as the projection 
method and similar techniques. This must be a final goal towards a real 
understanding of the true many-body dynamics. At this point, it is worth- 
while to underline that the above discussion is a heuristic scheme only but 
not a straightforward recipe. The specific method of introducing the IGFs 
depends on the form of operators An, the type of the Hamiltonian, and con- 
ditions of the problem. The irreducible parts in higher-order equations and 



connection with Mori formalism was considered by Tserkovnikov |51]. The 
incorporation of irreducible parts in higher-order equations and connection 
with the moment expansion was studied in ref. ( see Appendix B ). 
Here a sketchy form of the IGF method is presented. The aim to introduce 
the general scheme and to lay the groundwork for generalizations and spe- 
cific applications is expounded in the next Sections. We demonstrate below 
that the IGF method is a powerful tool for describing the quasi-particle 
excitation spectra, allowing a deeper understanding of elastic and inelastic 
quasi-particle scattering effects and corresponding aspects of damping and 
finite lifetimes. In the present context, it provides a clear link between the 
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equation-of-motion approach and the diagrammatic methods due to deriva- 
tion of the Dyson equation (37). Moreover, due to the fact that it allows the 
approximate treatment of the self-energy effects on a final stage, it yields a 
systematic way of the construction of approximate solutions. 
It is necessary to emphasize that there is an intimate connection between an 
adequate introduction of mean fields and internal symmetries of the Hamil- 
tonian. To test these ideas further, in the following Sections, we analyze 
the mean field and generalized mean field concepts for various many-body 
systems on a lattice. 

4 Many-Particle Interacting Systems on a Lattice 
4.1 Spin Systems on a Lattice 

There exists a big variety of magnetic materials. The group of magnetic 
insulators is of a special importance. For the group of systems considered 
in this Section, the physical picture can be represented by a model in which 
the localized magnetic moments originating from ions with incomplete shells 
interact through a short-range interaction. Individual spin moments form 
a regular lattice. The first model of a lattice spin system was constructed 
to describe a linear chain of projected electron spins with nearest-neighbor 
coupling. This was the famous Lenz-Izing model which was thought to yield 
a more sophisticated description of ferromagnetism than the Weiss uniform 
molecular field picture. However, in this model, only one spin component is 
significant. As a result, the system has no collective dynamics. The quantum 
states that are eigenstates of the relevant spin components are stationary 
states. The collective dynamics of magnetic systems is of great importance 
since it is related to the study of low-lying excitations and their interactions. 
This is the main aim of the present consideration. Although the Izing model 
was an intuitively right step forward from the uniform Weiss molecular field 
picture, the physical meaning of the model coupling constant remained com- 
pletely unclear. The concept of the exchange coupling of spins of two or more 
nonsinglet atoms appeared as a result of the Heitler-London consideration 
of chemical bond. This theory and the Dirac analysis of the singlet-triplet 
splitting in the helium spectrum stimulated Heisenberg to make a next es- 
sential step. Heisenberg suggested that the exchange interaction could be 
the relevant mechanism responsible for ferromagnetism. 

4.1.1 Heisenberg Ferromagnet 

The Heisenberg model of a system of spins on various lattices ( which was 
actually written down explicitly by van Vleck ) is termed the Heisenberg fer- 
romagnet and establishes the origin of the coupling constant as the exchange 
energy. The Heisenberg ferromagnet in a magnetic field H is described by 
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the Hamiltonian 

(39) H = -Y,J{i- j)S,S, - giisH J2 Sf 

ij i 

The couphng coefficient J{i — j) is the measure of the exchange interaction 
between spins at the lattice sites i and j and is defined usually to have 
the property J(i - j = 0) = 0. This constraint means that only the inter- 
exchange interactions are taken into account. However, in some complicated 
magnetic salts, it is necessary to consider an "effective" intra-site (see fS^ ) 
interaction (Hund-rule-type terms). The coupling, in principle, can be of a 
more general type (non-Heisenberg terms). These aspects of construction of 
a more general Hamiltonian are very interesting, but we do not pause here 
to give the details. 

For crystal lattices in which every ion is at the centre of symmetry, the 
exchange parameter has the property 

J{i-j) = J{j-i) 
We can rewrite then the Hamiltonian (|39| ) as 

(40) H = -J2Ji^- MSfS^ + StSr) 

ij 

Here = i: iS^ are the raising and lowering spin angular momentum 
operators. The complete set of spin commutation relations is 

[St,Sj]^ = 2St5i,; [St,Sr]^ = 2S{S + 1) - 2(5f)2; 
[Sf, = ±Sf6,,; St = S{S + 1) - {S[f - SrSt; 

(5+)25+i = 0, (5-)2^+i = 

We omit the term of interaction of the spin with an external magnetic field 
for the brevity of notation. The statistical mechanical problem involving 
this Hamiltonian was not exactly solved, but many approximate solutions 
were obtained. 

To proceed further, it is important to note that for the isotropic Heisenberg 
model, the total z-component of spin S^^^ = J2i is a constant of motion, 
i.e. 

[H,SU = 

There are cases when the total spin is not a constant of motion, as, for 
instance, for the Heisenberg model with the dipole terms added. 
Let us define the eigenstate iV'o > so that S^ltpo >= for all lattice sites 
Ri. It is clear that {ipo > is a state in which all the spins are fully aligned 
and for which Sf\ipo >= S\ipo >. We also have 

i 
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, where the reciprocal vectors k are defined by cychc boundary conditions. 
Then we obtain 



Here N is the total number of ions in the crystal. So, for the isotropic 
Heisenberg fcrromagnet, the ground state \ipo > has an energy —NS^Jq. 
The state [ipo > corresponds to a total spin NS. 

Let us consider now the first excited state. This state can be constructed 
by creating one unit of spin deviation in the system. As a result, the total 
spin is NS — 1. The state 

is an eigenstate of R which corresponds to a single magnon of the energy 
(41) iJi'^\k)=2S{J^-J^) 

Note that the role of translational symmetry, i.e. the regular lattice of spins, 
is essential, since the state iV'jt > is constructed from the fully aligned state 
by decreasing the spin at each site and summing over all spins with the 
phase factor e^^^i . It is easy to verify that 

<'^k\SloMk >=NS-\ 



The above consideration was possible because we knew the exact ground 
state of the Hamiltonian . There are many models where this is not the 
case. For example, we do not know the exact ground state of a Heisenberg 
ferromagnet with dipolar forces and the ground state of the Heisenberg 
antiferr omagnet . 

4.1.2 Heisenberg Antiferromagnet 

We now discuss the Heisenberg model of the antiferromagnet which is more 
complicated to analyse. The fundamental problem here is that the exact 
ground state is unknown. We consider, for simplicity, a two-sublattice struc- 
ture in which nearest neighbour ions on opposite sublattices interact through 
the Heisenberg exchange. For a system of ions on two sublattices, the Hamil- 
tonian is 

(42) = J E SmSjn+5 + E ^nSn+S 

m,5 n,S 

Here the notation m = Rm means the position vectors of ions on one 
sublattice (a) and n for the ions on the other (5). Nearest neighbor ions 
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on different sublattices are a distance \6\ apart. ( The anisotropy field 
IJ-HAiJ2m ~ Sn'S'n)) whicli is not Written down explicitly, is taken to 
be parallel to the z axis. ) The simplest crystal structures that can be 
constructed from two interpenetrating identical sublattices are the body- 
centered and simple cubic. 

The exact ground state of this Hamiltonian is not known. One can use the 
approximation of taking the ground state to be a classical ground state, 
usually called the Neel state, in which the spins of the ions on each sub- 
lattice are oppositely aligned along the z axis. However, this state is not 
even an eigenstate of the Hamiltonian (^). Let us remark that the total 
z-component of the spin commutes with the Hamiltonian (^). It would be 
instructive to consider here the construction of a spin wave theory for the 
low-lying excitations of the Heisenberg antiferromagnet in a sketchy form to 
clarify the foregoing. 

To demonstrate the specifics of Heisenberg antiferromagnet more explicitly, 
it is convenient to rotate the axes of one sublattice through tt about the 
X-axis. This transformation preserves the spin operator commutation rela- 
tions and therefore is canonical. Let us perform the transformation on the 
Rn, or 6-sublattice 

qz qz . c<± cT 

The operators S"^ and S!^ commute, because they refer to different sublat- 
tices. 

The transformation to the momentum representation is modified in compar- 
ison with the ferromagnet case 

q q 
Here q is the reciprocal lattice vectors for one sublattice, each sublattice 



containing N ions. After these transformations, the Hamiltonian (42) can 
be rewritten as 

(43) H = ^Y1 2^^5[(5-5+ + S-S+) + j,{S+S+ + 5-5")] 

In (^3D, 7q is defined as zjg = J2m=n.n. ^xp{iqRm), and z is the number of 
nearest neighbors; the constant terms and the products of four operators are 
omitted. Thus the Hamiltonian of the Heisenberg antifferomagnet is more 
complicated than that for the ferromagnet. Because it contains two types 
of spin operators that are coupled together, the diagonalization of ( ^3| ) has 
its own specificity. 

To diagonalize (|4^), let us make a linear transformation to new operators ( 
Bogoliubov transformation ) 

(44) 5+ = UqUg + Vqh\\ S~ = Uqh\ + Vqttq 



19 



with 

The transformation coefficients and Vk are purely real. To preserve the 
commutation rules for the spin operators 

[S+,S^,]=2SN6k,k' 

, they should satisfy u'^{k) — v'^{k) = 2SN. The transformations from the 
operators {S^,S~) to the operators (ag,6j) give 

(ajag + blbg)[{u^{q) + v^{q)) + 2uqVg-iq] 
+ {aghg + a\b\)[{u^{q) + v^{q))-iq + 2ugVg\ 
(45) +2ugVg-ig + 2v'^{q) 



We represented Hamiltonian ( ^3] ) as a form quadratic in the Bose opera- 
tors (ag,6j) . We shall now consider the problem of diagonalization of this 



form|^6[. To diagonalize (43), we should require that 

2uqVq + {u^{q)+v'^{q))-1g = ^ 

Then we obtain 

(46) 2u^{q) = 25ivii±^; 2v'^{q) = 2SN 



Kg 



Here the following notation was introduced: Kg = ^(1 — 7^) and 2uqVq 
—2SNjq/ Kg After the transformation (^^, we get, instead of (|4^) , 



(47) H = J24^f"'\k){alaq + blbg 

k 

with 



(48) <^"'^(^) = 2.J5v/l-7| 



Expression ( ^7D contains two terms, each with the same energy spectrum. 
Thus, there are two degenerate spin wave modes, because there can be 
two kinds of precession of the spin about the anisotropy direction. The 
degeneracy is lifted by the application of an external magnetic field in the 
z direction, because in this case the two sublattices become nonequivalent. 
These results should be kept in mind when discussing the quasi-particle 
many-body dynamics of the spin lattice models. 
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4.2 Correlated Electrons on a Lattice 

The importance of intra-atomic correlation effects in determining the mag- 
netic properties of transition metals and their compounds and oxides was 
recognized many years ago. The essential basis of studies of metallic mag- 
netism, namely, that the dominant physical mechanism responsible for the 
observed magnetic properties of the transition metals and their compounds 
and alloys is the strong intra-atomic correlation in an otherwise tight-binding 
picture, is generally accepted as being most suitable. The problem of the 
adequate description of strongly correlated electron systems on a lattice 
was studied intensively during the last decade, especially in the context of 
metallic magnetism, heavy fermions, and high- Tc superconductivity 0. The 
understanding of the true nature of electronic states and their quasi-particle 
dynamics is one of the central topics of the current experimental and the- 
oretical efforts in the field. The source of spin magnetism in solids is, of 
course, the Pauli exclusion principle as manifested in the exchange interac- 
tion and higher order mechanism. Of particular interest is the fact that the 
Hartree-Fock or mean field theory, i.e. the theory including exchange but 
not correlation effects, invariably overestimates the tendency to magnetism. 
This fact obviously complicated the already complicated problem of mag- 
netism in a metal with the d band electrons which, as was mentioned above, 
are really neither "local" nor "itinerant" in a full sense. 
The strongly correlated electron systems are systems in which electron cor- 
relations dominate. The theoretical studies of strongly correlated systems 
had as a consequence the formulation of two model Hamiltonians which 
play a central role in our attempts to get an insight into this complicated 
problem. These are the Anderson single-impurity model (SIAM) [^] and 
Hubbard model |^^. It was only relatively recently recognized that both 
the models have a very complicated many-body dynamics, and their "sim- 
plicity" manifests itself in the dynamics of two-particle scattering, as was 
shown via elegant Bethe-anzatz solutions. 

4.2.1 Hubbard Model 

The model Hamiltonian usually referred to as the Hubbard Hamiltonian 
(49) ^ = Y1 + X! '^i<^'^i-^ 

ija id 

includes the intra-atomic Coulomb repulsion U and the one-electron hopping 
energy tij . The electron correlation forces electrons to localize in the atomic 
orbitals which are modelled here by a complete and orthogonal set of the 
Wannier wave functions [(f){f — Rj)]. On the other hand, the kinetic energy 
is reduced when electrons are delocalized. The main difficulty in solving 
the Hubbard model correctly is the necessity of taking into account both 



21 



these effects simultaneously. Thus, the Hamiltonian (^) is specified by two 
parameters: U and the effective electron bandwidth 



,2x1/2 



A = iN-'J2\U 

The band energy of Bloch electrons e{k) is defined as follows 

= N-^ ^i^) exp[ik{Ri - Rj], 

k 

where is the number of lattice sites. It is convenient to count the energy 
from the center of gravity of the band, i.e. ta = to = Sfc^(^) = ( 
sometimes it is useful to retain to explicitly ) . 

This conceptually simple model is mathematically very complicated. The 
effective electron bandwidth A and Coulomb intra-site integral U determine 
different regimes in 3 dimensions depending on the parameter 7 = A/U. In 
addition, the Pauli exclusion principle that does not allow two electrons of 
common spin to be at the same site, i.e. nf^ = riia-, plays a crucial role, and it 
should be taking into account properly while making any approximations. It 
is usually rather a difficult task to find an interpolating solution for dynamic 
properties of the Hubbard model for various mean particle densities. To 
solve this problem with a reasonably accuracy and to describe correctly an 
interpolated solution from the "band" limit (7 ^ 1) to the "atomic" limit 
(7 0), one needs a more sophisticated approach than usual procedures 
developed for description of the interacting electron gas problem [^]. We 
have evidently to improve the early Hubbard theory taking into account of 
variety of possible regimes for the model depending on the electron density, 
temperature, and values of 7. The single-electron GF 

(50) Gijaiuj) =« ai^\a]^ »= A^"^ ^ w) exp[-?A?(^i - ^j)]. 



calculated by Hubbard [54|, [55|, has the characteristic two-pole functional 
structure 

(51) G^{k,Lo) = [F^{uj)-e{k)]-^ 
where 

(52) F„ (w) = — — ^— + _ a 

Here = n , = 1 — n; Ej^ = U, -E_ = 0, and A is a certain function 
which depends on parameters of the Hamiltonian. In this approximation, 
Hubbard took account of the scattering effect of electrons with spins a by 
electrons with spin —a which are frozen as well as the "resonance broaden- 
ing" effect due to the motion of the electrons with spin —a. The "Hubbard 
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Ill" decoupling procedure suffered of serious limitations. However, in spite 
of the limitations, this solution gave the first clue to the qualitative un- 
derstanding of the property of narrow-band system like the metal-insulator 
transition. 

If A is small (A — > 0), then expression (^2|) takes the form: 



uj — -E_ — nlg-A u) — — n_^X ' 
which corresponds to two shifted subbands with the gap 

(Ji - W2 = {E+ - E.) + (nl^ - nt^)\ = U + X2nt^. 
If A is very big, then we obtain 

A 1 



[(a; - E.)nZ^ + {to - E+)nt^]X oo - {nt,E+ - nZ^E. 



The latter solution corresponds to a single band centered at the energy 
(jO ~ n'tfjU. Thus, this solution explains qualitatively the appearance of a 
gap in the density of states when the value of the intra-atomic correlation 
exceeds a certain critical value, as it was first conjectured by N. Mott. 
The two-pole functional structure of the single-particle GF is easy to under- 
stand within the formalism that describes the motion of electrons in binary 
alloys (5^,|6^]. If one introduces the two types of the scattering potentials 
t± ^ {uj — E-i-)~^, then the two kinds of the t-matrix and T_ appears 
which satisfy the following system of equations: 



where is the bare propagator between the sites with energies E±. The 
solution of this system is of the following form 

t± + t±Git± 



(53) 



(1 - t+G\+)(l - t^G^__) - G%Gl_t+t^ 

t^' + Gj 

(t;^ - Gl+){tZ^ - _) - G^_+Gl_ ' 



Thus, by comparing this functional two-pole structure and the "Hubbard 
HI" solution 

= UJ - F^{uj) 
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, it is possible to identify the "scattering corrections" and "resonance broad- 
ening corrections" in the following way: 



uj{(jj — U) — {uj — Un^a)A(jioj) 



A^uj) = Y^uj) + Y^^uj) - Yl^U - to) 



If we put A„{uj) = 0, we immediately obtain the "Hubbard I" solution ||5J1 



(54) Gp)(fc,a;)« -^^^ + ^ " 



UJ — U — e{k)n^fj uj — e{k){l — n^a) 

Despite that this solution is exact in the atomic limit ( tij = 0), the "Hub- 
bard I" solution has many serious drawbacks. The corresponding spectral 
function consists of two (5- function peaks. The "Hubbard HI" solution in- 
cludes several corrections, including scattering corrections which broadens 
the peaks and shift them when U is changed. 

The "alloy analogy" approximation corresponds to A(j[ijo) ~ l^(i^). An in- 
teresting analysis of the "Hubbard HI" solution was performed in paper [|6C|]. 
The Hubbard sub-band structure was obtained in an analytic form in the 
"Hubbard HI" approximation, using the Lorentzian form for the density of 
states for non-interacting electrons. This resulted in an analytical form for 
the self-energy and the density of states for interacting electrons. Note that 
the "Hubbard HI" self-energy operator is local, i.e. does not depend 

on the quasi-momentum. Another drawback of this solution is a very incon- 
venient functional representation of elastic and inelastic scattering processes. 
The conceptually new approach to the theory of very strong but finite elec- 



tron correlation for the Hubbard model was proposed by Roth [9C]. She 
clarified microscopically the origination of the two-pole solution of the single- 
particle GF in the strongly correlated limit 



(55) G(«)(fe, 



UJ] 



ri-rr 1 — n_ 



UJ - U - e{k)n-a - Wfc_<x(l - n_<^) w - e{k){l - n_„) - n_„Wk-a 

We see that, in addition to a band narrowing effect, there is an energy shift 
Wk-a given by 

ncr(l - n„)Wka = X! *u < - 'f^i-cr - n-j-a) > - 51 exp[iA;(j - i)] 

ij ij 

(56) {nl- < Ui^rija > + < aj_^a|^aj<:rat_^ > + < a^j^^a^j^aiaai^^ >) 
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This energy shift corrects the situation with the "Hubbard I" spectral func- 
tion and recovers, in principle, the possibility of describing the ferromagnetic 
solution. Thus, the Roth solution gives an improved version of "Hubbard 
I" two-pole solution and includes the band shift, that is most important 
in the case of a nearly-half-filled band. It is worth noting that this result 
was a very unusual fact from the point of view of the standard Fermi-liquid 
approach, showing that the naive one-electron approximation of band struc- 
ture calculations is not valid for the description of electron correlations of 
lattice fermions. 

It is this feature - the strong modification of single-particle states by many- 
body correlation effects - whose importance we wish to emphasize here. 
Various attempts were made to describe the properties of the Hubbard model 
in both the strong and weak coupling regimes and to find a better solution 
{e.g. 1 56] - ||5^ ). Different schemes of decoupling of the equations of 



motion for the GFs analysed and compared in paper ||59[], when calculating 
the electron contribution to the cohesive energy in a narrow band system. 
These calculations showed importance of the correlation effects and the right 
scheme of approximation. 

Thus, a sophisticated many-body technique is to be used for calculating the 
excitation spectra and other characteristics at finite temperatures. We shall 
show here following papers [43|,[23| that the IGF method permits us to im- 



prove substantially both the solutions, Hubbard and Roth, by defining the 
correct Generalized Mean Fields for the Hubbard model. 

4.2.2 Single Impurity Anderson Model (SIAM) 

The Hamiltonian of SIAM can be written in the form [p3l 



(57) H = Y^ ^kcl^Cka + ^OafLfoa + f//2 ^ riQ^no^^ + 

a a 



ka 



where c^^ and /J^ are, respectively, the creation operators for conduction 
and localized electrons; is the conduction electron energy, £"00- is the lo- 
calized electron energy level, and U is the intra-atomic Coulomb interaction 
at the impurity site; Vk represents the s — {d)f hybridization interaction 



term and was written in paper 1 53 1 in the following form 

1 

The hybridization matrix element is 



(58) Vk = -^J2 Vf{Rj)exp{ikRj 
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The use of Hartree-Fock term here is notable, since it justifies the initial 



treatment of SI AM in[53| entirely in the H-F approximation. A number 
of approaches for SIAM and other correlated electronic systems was pro- 
posed, aimed at answering the Anderson question: "...whether a real many- 
body theory would give answer radically different from the Hartree-Fock 
results?" H. 

Our goal is to propose a new combined many-body approach for the descrip- 
tion of many-body quasi-particle dynamics of SIAM at finite temperatures. 
The interplay and competition of the kinetic energy (ek), potential energy 
(U), and hybridization (V) substantially infiuence the electronic spectrum. 
The renormalized electron energies are temperature-dependent, and elec- 
tronic states have finite lifetimes. These effects are described most suitable 
by the Green functions method. The purpose of the present approach is to 
find the electronic quasi-particle spectrum renormalized by interactions (U- 
and V-terms) in a wide range of temperatures and model parameters and 
to calculate explicitly the damping of the electronic states. 

4.2.3 Periodic Anderson Model (PAM) 

Let us now consider a lattice generalization of SIAM, the so-called periodic 
Anderson model (PAM). The basic assumption of the periodic impurity 
Anderson model is the presence of two well-defined subsystems, i.e. the 
Fermi sea of nearly free conduction electrons and the localized impurity 
orbitals embedded into the continuum of conduction electron states ( in 
rare-earth compounds, for instance, the continuum is actually a mixture of 
s, p, and d states, and the localized orbitals are / states). The simplest form 
of PAM 

(59) H = Y^ ^kc\,cka + E ^o4/ia + ^7/2 Ui^m., + 

ka ia ia 

V 



Y{exp{ikRi)c\^fi^ + exp{-ikRi)fl^Cka) 



^^ka 



assumes a one-electron energy level Eq, hybridization interaction V, and the 
Coulomb interaction U at each lattice site. Using the transformation 



^ ~ ^ e-xp{-^^Rjy]a'^ = -L exp{ikRj)cja 

j ^ 3 



the Hamiltonian ( |5^ ) can be rewritten in the Wannier representation: 

(60) H = Y tidaCjo + E Eoflfia + U/2 J2 niaUi-a + 

ia ia 

^E(4/- + /lc..) 



ija 
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If one retains the A;-dependence of the hybridization matrix element 14 in 



(60) , the last term in the PAM Hamiltonian describing the hybridization 
interaction between the localized impurity states and extended conduction 
states and containing the essence of a specificicity of the Anderson model, 
is as follows 

^ + flc,,y, V,, = ^ ^ Vkexp[ik[R, - Ri)] 

ijcr k 

The on-site hybridization Va is equal to zero for symmetry reasons. A 
detailed analysis of the hybridization problem from a general point of view 
and in the context of PAM was made in paper 1 61]. The Hamiltonian of 
PAM in the Bloch representation takes the form 

(61) H = J2 efecLcfc. + E Ekflfka + U/2 ^ m^m.^ + 

ka icr icr 

^kicl^fka + fLcka) 

ka 

Note that as compared to the SIAM, the PAM has its own specific features. 
This can lead to peculiar magnetic properties for concentrated rare-earth 
systems where the criterion for magnetic ordering depends on the competi- 
tion between indirect RKKY-type interaction p^] ( not included into SIAM ) 
and the Kondo-type singlet-site screening ( contained in SIAM ). The inclu- 
sion of inter-impurity correlations makes the problem more difficult. Since 
these inter-impurity effects play an essential role in physical behaviour of 
real systems ||6^ , ||63[] , it is instructive to consider the two-impurity Anderson 
model (TIAM) too. 

4.2.4 Tvi^o-Impurity Anderson Model ( TIAM ) 



The two- impurity Anderson model was considered by Alexander and Anderson 1 64] 
They put forward a theory which introduces the impurity-impurity interac- 
tion within a game of parameters. The Hamiltonian of TIAM reads 

(62) H = Y,iijc\^Cj,+ EQiflU, + U/2 ^ n.^n^., + 

ija i=l,2a i=l,2a 

J2(^k^clfia + V,kflc^.) + J2(^12flf2a + V2lflfla) 
ia a 

where Eqi are the position energies of localized states ( for simplicity, we con- 
sider identical impurities and s-type (i.e. non-degenerate) orbitals: £^01 = 
Eq2 = Eq. Let us recall that the hybridization matrix element Vik was de- 
fined in ([5^). As for the TIAM, the situation with the right definition of the 
parameters V12 and Vik is not very clear. The definition of V12 inj^] is the 
following 

V12 = Vl = J 4{r)HfM^dr 
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(now Hf without "H-F" mark). The essentiahy local character of the Hamil- 
tonian Hf clearly shows that V12 describes the direct coupling between near- 
est neighboring sites ( for a detailed discussion see[29| where the hierarchy 
of the Anderson models was discussed too ). 



5 Effective and Generalized Mean Fields 
5.1 Molecular Field Approximation 

The most common technique for studying the subject of interacting many- 
particle systems is to use the mean field theory. This approximation is 



especially popular in the theory of magnetism |65]. Nevertheless, it was 
pointed |Q that 



"the Weiss molecular field theory plays an enigmatic role in sta- 
tistical mechanics of magnetism" . 

To calculate the susceptibility and other characteristic functions of a system 
of localized magnetic moments, with a given interaction Hamiltonian, the 
approximation, termed the " molecular field approximation" was used widely. 
However, it is not an easy task to give the formal unified definition what 
the mean field is. In a sense, the mean field is the umbrella term for a 
variety of theoretical methods of reducing the many-particle problem to the 
single-particle one. Mean field theory, that approximates the behaviour of 
a system by ignoring the effect of fiuctuations and those spin correlations 
which dominate the collective properties of the ferromagnet usually provides 
a starting and estimating point only, for studying phase transitions. The 
mean field theories miss important features of the dynamics of a system. 
The main intention of the mean field theories, starting from the works of van 
der Waals and P.Weiss, is to take into account the cooperative behaviour 
of a large number of particles. It is well known that earlier theories of 
phase transitions based on the ideas of van der Waals and Weiss lead to 
predictions which are qualitatively at variance with results of measurements 
near the critical point. Other variants of simplified mean field theories such 
as the Hartree-Fock theory for electrons in atoms, etc lead to discrepancies 
of various kinds too. It is therefore natural to analyze the reasons for such 
drawbacks of earlier variants of the mean field theories. 



5.2 Effective Field Theories 

A number of effective field theories which are improved versions of the 
"molecular field approximation" were proposed. It is the purpose of this 
study to stress a specificity of strongly correlated many-particle systems on 
a lattice contrary to continuum (uniform) systems. Although many impor- 
tant questions remain still unresolved, a vision of useful synthesis begins 
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to emerge. As a workable eye-guide , the set of mean field theories ( most 
probably incomplete ) is shown in Table 1. The meaning of many these 
entries and terms will become clearer in the forthcoming discussion and will 
put them in a clearer perspective. My main purpose is to elucidate ( at 
least in the mathematical structure ) and to give plausible arguments for 
the tendency, which expounded in Table 1. This tendency shows the fol- 
lowing. The earlier concepts of molecular field were described in terms of a 
functional of mean magnetic moments (in magnetic terminology ) or mean 
particle densities ( Hartree-Fock field ) . The corresponding mean-field func- 
tional F[< n >, < 5^ >] describes the uniform mean field. 
Actually, the Weiss model was not based on discrete "spins" as is well known, 
but the uniformity of the mean internal field was the most essential feature 
of the model. In the modern language, one should assume that the interac- 
tion between atomic spins Si and its neighbors is equivalent to a mean ( or 
molecular ) field. Mi = Xo[^i'^^*'' + ^i"^'^'^] that the molecular field h["^'^^ 
is of the form /i^"^-^) = J^i J{Rji) < Si > (above Tc ). Here /i*^^* is an applied 
conjugate field, xo is the response function, and J{Rji) is an interaction. In 
other words, the mean field approximation reduces the many-particle prob- 
lem to a single-site problem in which a magnetic moment at any site can 
be either parallel or antiparallel to the total magnetic field composed of the 
applied field and the molecular field. The average interaction of i neigh- 
bors was taken into account only, and the fluctuations were neglected. One 
particular example, where the mean fleld theory works relatively well is the 
homogeneous structural phase transitions; in this case the fluctuations are 
conflned in phase space. 

The next important step was made by L. Neel |6?]. He conjectured that 
the Weiss internal fleld might be either positive or negative in sign. In the 
latter case, he showed that below a critical temperature ( Neel temperature 
) an ordered arrangement of equal numbers of oppositely directed atomic 
moments could be energetically favorable. This new magnetic structure was 
termed antiferromagnetism. It was conjectured that the two-sublattice Neel 
( classical ) ground state is formed by local staggered internal mean flelds. 
There is a number of the "correlated effective fleld" theories, that tend to 
repair the limitations of simplifled mean fleld theories. The remarkable and 
ingenious one is the Onsager "reaction fleld approximation" 1 68]. He sug- 
gested that the part of the molecular fleld on a given dipole moment which 
comes from the reaction of neighboring molecules to the instantaneous ori- 
entation of the moment should not be included into the effective orienting 
fleld. This "reaction fleld" simply follows the motion of the moment and 
thus does not favor one orientation over another. The meaning of the mean 
fleld approximation for the spin glass problem is very interesting but spe- 
ciflc, and we will not discuss it here. A single-site molecular-fleld model 
for randomly dilute ferro- and antiferromagnets in the framework of the 
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double-time thermal GFs was presented in paper [|69[. 



5.3 Generalized Mean Fields 



It was shown [39|, |^^, |7C] that mean- field approximations, for example the 
molecular field approximation for a spin system, the Hartree-Fock approxi- 
mation and the BCS-Bogoliubov approximation for an electron system are 
universally formulated by the Peierls-Bogoliubov-Feynman (PBF) inequal- 
ity: 

Here F is the free energy, and is a " trial" or a " mean field" approximat- 
ing Hamiltonian. This inequality gives the upper bound of the free energy of 
a many-body system. It is important to emphasize that the BCS-Bogoliubov 
theory of superconductivity ||l^,|7l| was formulated on the basis of a trial 
Hamiltonian which consists of a quadratic form of creation and annihila- 
tion operators, including "anomalous" ( off-diagonal ) averages |10|. The 
functional of the mean field ( for the superconducting single-band Hubbard 
model ) is of the following form [71|: 



(64) K = U 



. - < 4-a4a > - < Aa^i^ > 



The "anomalous" off-diagonal terms fix the relevant BCS-Bogoliubov vac- 
uum and select the appropriate set of solutions. 

Another remark about the BCS-Bogolubov mean-field approach is instruc- 
tive. Speaking in physical terms, this theory involves a condensation cor- 
rectly, in spite that such a condensation cannot be obtained by an expansion 
in the effective interaction between electrons. Other mean field theories, e.g. 
the Weiss molecular field theory and the van der Waals theory of the liquid- 
gas transition are much less reliable. The reason why a mean- field theory of 
the superconductivity in the BCS-Bogoliubov form is successful would ap- 
pear to be that the main correlations in metal are governed by the extreme 
degeneracy of the electron gas. The correlations due to the pair condensa- 
tion, although they have dramatic effects, are weak ( at least in the ordinary 
superconductors ) in comparison with the typical electron energies, and may 
be treated in an average way with a reasonable accuracy. All above remarks 
have relevance to ordinary low-temperature superconductors. In high-Tc 
superconductors, the corresponding degeneracy temperature is much lower, 
and the transition temperatures are much higher. In addition, the relevant 
interaction responsible for the pairing and its strength are unknown. From 
this point of view, the high-Tc systems are more complicated. It should be 
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clarified what governs the scale of temperatures, i.e. critical temperature, 
degeneracy temperature, interaction strength or their complex combination, 
etc. In this way a useful insight into this extremely complicated problem 
would be gained. 

Generalization of the molecular field approximation on the basis of the PBF 
inequality is possible when we know a particular solution of the model {e.g., 
for one-dimensional Ising model we know the exact solution in the field). 
One can use this solution to get a better approximation than the mean field 
theory. There are some other methods of improvement of the molecular field 
theory |72|, [ff^]. Unfortunately, these approaches are nonsystematic. 



From the point of view of quantum many-body theory, the problem of ad- 
equate introduction of mean fields for system of many interacting particles 
can be most consistently investigated in the framework of the IGF method. 
A correct calculation of the quasi-particle spectra and their damping, partic- 



ularly, for systems with a complicated spectrum and strong interaction |24] 
reveals, as it will be shown below, that the generalized mean fields can have 
very complicated structure which cannot be described by a functional of the 
mean-particle density. 

To illustrate the actual distinction of description of the generalized mean 
field in the equation-of-motion method for the double-time Green functions. 



let us compare the two approaches, namely, that of Tyablikov |46] and of 



Callen We shall consider the Green function << S^\S » for the 

isotropic Heisenberg model 



(65) ^ = 

ij 

The equation of motion (|T^ ) for the spin Green function is of the form 

(66) uj « StlSj- »^= 
2 < 5" > 6ij + ^ J(i - g) « St SI - StS!\Sj »^ 

9 

The Tyablikov decoupling expresses the second-order GF in terms of the 
first (initial) GF: 

(67) « StSl\Sj »=< S^ >« St\Sj- » 

This approximation is an RPA-type; it does not lead to the damping of spin 
wave excitations (c/. (|4l|) ) 



E{q) = Y.J{i-g)<S'> exp[i{Ri - Rg)q\ = 2 < 5^ > ( Jq - J,) 

a 

The reason for this is rather transparent. This decoupling does not take into 
account the inelastic magnon-magnon scattering processes. In a sense, the 
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Tyablikov approximation consists of approximating the commutation rela- 
tions of spin operators to the extent of replacing the commutation relation 
[St, S.]. = 2SfS,, by [St, S-]_ = 2<S^> 5,, . 



Callen |74| has proposed an improved decoupling approximation in the 



method of Tyablikov in the following form: 

(69) « SlSj\B »^< 5^ >« S+\B » -a < Sg S^ >« 5+|B » 

Here < a < 1. To clarify this point, it should be reminded that for spin 
1/2 ( the procedure was generalized by Callen to an arbitrary spin), the spin 
operator S^ can be written as Sg = S — Sg S^ or 5| = ^('S'^'S'" — S~St)- 
It is easy to show that 

1 — (y. 1 ~|~ (y. 

Sg = aS -\ ^9 ^9 

The operator SgS^ represents the deviation of < 5^ > from S. In the 
low-temperature region, this deviation is small, and a ~ 1. Similarly, the 
operator ^{St S~ — S" St) represents the deviation of < 5^ > from 0. Thus, 
when < S^ > approaches to zero, one can expect that a ^ 0. Thus, in this 
way, it is possible to obtain a correction to the Tyablikov decoupling with 
either a positive or negative sign, or no correction at all, or any intermediate 
value, depending on the choice of a. The above Callen arguments are not 
rigorous , for, although the difference in the operators S'^S' and S~ S^ is 
small if < S^ 0, each operator makes a contribution of the order of 5, and 
it is each operator which is treated approximately, not the difference. There 
are some other drawbacks of the Callen decoupling scheme. Nevertheless, 
the Callen decoupling was the first conceptual attempt to introduce the 
interpolation decoupling procedure. Let us note that the choice of a = 
over the entire temperature range is just the Tyablikov decoupling (67). 
The energy spectrum for the Callen decoupling is given by 

(70) E(5) = 2 < 5^ > ((Jo - J,) + E['^(^) - - QW(.E{k))) 

k 

Here N{E{k)) is the Bose distribution function N{E{k)) = [exp{E{k)l3) - 
This is an implicit equation for N{E{k)), involving the unknown quan- 
tity < 5^ > . For the latter an additional equation is given ||7^. Thus, both 
these equations constitute a set of coupled equations which must be solved 
self-consistently for < S^ >. 

This formulation of the Callen decoupling scheme displays explicitly the ten- 
dency of the improved description of the mean field. In a sense, it is possible 
to say that the Callen work dates really the idea of the generalized mean 
field within the equation-of-motion method for double-time GFs, however. 



in a semi- intuitive form. The next essential steps were made by Plakida |47] 



for the Heisenberg ferromagnet and by Kuzemsky 143], |23] for the Hubbard 
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model. As was mentioned above, the correct definition of Generalized Mean 
Fields depends on the condition of the problem, the strength of interaction, 
the choice of relevant operators, and on the symmetry requirements. 



5.4 Symmetry Broken Solutions 

In many-body interacting systems, the symmetry is important in classi- 
fying different phases and in understanding the phase transitions between 
them According to Bogoliubov |7|]( cf. refs. 0, in each 



condensed phase, in addition to the normal process, there is an anomalous 
process (or processes) which can take place because of the long-range in- 
ternal field, with a corresponding propagator. Additionally, the Goldstone 
theorem[^] states that, in a system in which a continuous symmetry is bro- 
ken (i.e. a system such that the ground state is not invariant under the 
operations of a continuous unitary group whose generators commute with 
the Hamiltonian ), there exists a collective mode with frequency vanishing, 
as the momentum goes to zero. For many-particle systems on a lattice, 
this statement needs a proper adaptation. In the above form, the Gold- 
stone theorem is true only if the condensed and normal phases have the 
same translational properties. When translational symmetry is also broken, 
the Goldstone mode appears at a zero frequency but at nonzero momen- 
tum, e.g., a crystal and a helical spin-density- wave (SDW) ordering (see for 
discussionQ-H). 

The anomalous propagators for an interacting many-fermion system corre- 
sponding to the ferromagnetic (FM) , antiferromagnetic (AFM), and super- 
conducting (SC) long-range ordering are given by 

(71) FM : Gfm ~« flfca; 4-a » 

AFM : Gafm ~<< «fc+Qf7;4+Q'a' >> 

SC : Gsc ~<< ctfco-; fl-fe-o- >> 



In the SDW case, a particle picks up a momentum Q — Q' from scattering 
against the periodic structure of the spiral ( nonuniform ) internal field, 
and has its spin changed from a to a' by the spin-aligning character of the 
internal field. The Long-Range-Order (LRO) parameters are: 

(72) FM ■.m = l/Nj2< Aa^^k-. > 

k<T 

AFM :Mq = J2< aia^k+Q-a > 

ka 

SC:A = J2< ali^^al^ > 
k 
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It is important to note that the long-range order parameters are functions 
of the internal field, which is itself a function of the order parameter. There 
is a more mathematical way of formulating this assertion. According to the 



paper |75], the notion "symmetry breaking" means that the state fails to 
have the symmetry that the Hamiltonian has. 

A true breaking of symmetry can arise only if there are infinitesimal " source 
fields". Indeed, for the rotationally and translationally invariant Hamilto- 
nian, suitable source terms should be added: 

(73) FM : EfiBH., E «L«fc-a 

fccr 

AFM : efLsH^al^ak+Q-a 

kQ 



SC : -fci^^fcT ^ 



fci"fcT ^ "fcT"-fci; 

k 

where e ^ is to be taken at the end of calculations. 

For example, broken symmetry solutions of the SDW type imply that the 
vector Q is a measure of the inhomogeneity or breaking of translational 
symmetry. The Hubbard model is a very interesting tool for analyzing the 
symmetry broken concept. It is possible to show that antiferromagnetic state 
and more complicated states ( e.g. ferrimagnetic) can be made eigenfunc- 
tions of the self-consistent field equations within an "extended" mean- field 
approach, assuming that the "anomalous" averages < al^ai-a > determine 
the behaviour of the system on the same footing as the "normal" density of 
quasi-particles < aj^aia >■ It is clear, however, that these "spin-flip" terms 
break the rotational symmetry of the Hubbard Hamiltonian. For the single- 
band Hubbard Hamiltonian, the averages < aj_^ai^cr >= because of the 
rotational symmetry of the Hubbard model. The inclusion of "anomalous" 
averages leads to the so-called " unresricted" H-F approximation (UHFA). 
This type of approximation was used sometimes also for the single-band 
Hubbard model for calculating the density of states. For this aim, the fol- 
lowing definition of UHFA 

(74) ni_„ai„ > ai„- < a\_^ai^ > aj_^ 

was used. Thus, in addition to the standard H-F term, the new so-called 
"spin-flip" terms are retained. This example clearly shows that the struc- 
ture of mean field follows from the specificity of the problem and should be 
defined in a proper way. So, one needs a properly defined effective Hamilto- 



nian -ffeflf- In paper |83| we thoroughly analyzed the proper definition of the 
irreducible GFs which includes the "spin-flip" terms for the case of itinerant 
antiferromagnetismp3] of correlated lattice fermions. For the single-orbital 



Hubbard model, the definition of the "irreducible" part should be modified 
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in the following way: 

{ir) t I t _ t it 

(75) 6pfl < Ug^a- > Gka- < flfc+pcrap+g-CT ><< aq-a\ala- »L0 

From this definition it follows that this way of introduction of the IGF 
broadens the initial algebra of operators and the initial set of the GFs. This 
means that the "actual" algebra of operators must include the spin-flip terms 
from the beginning, namely: (ojo-, a|^, njo-, a\^ai-o-). The corresponding 
initial GF will be of the form 

<< aia\a)j^ » « ai^|a]_^ >> \ 
<< ai-.a\a]^ » « ai-a\a]-a » / 

With this definition, one introduces the so-called anomalous (off-diagonal) 
GFs which fix the relevant vacuum and select the proper symmetry broken 
solutions. In fact, this approximation was investigated earlier by Kishore 
and Joshi |8^. They clearly pointed out that they assumed a system to be 
magnetised in the x direction instead of the conventional z axis. 
The problem of finding the ferromagnetic and antiferromagnetic "symme- 
try broken" solutions of the correlated lattice fermion models within IGF 
method was investigated in ref. |8^. A unified scheme for the construction 
of Generalized Mean Fields ( elastic scattering corrections ) and self-energy ( 
inelastic scattering ) in terms of the Dyson equation was generalized in order 
to include the " source fields" . The " symmetry broken" dynamic solutions of 
the Hubbard model which correspond to various types of itinerant antiferro- 
magnetism were discussed. This approach complements previous studies of 
microscopic theory of the Heisenberg antiferromagnet and clarifies the 
concepts of Neel sublattices for localized and itinerant antiferromagnetism 
and "spin-aligning fields" of correlated lattice fermions. 



6 Quasi-Particle Many Body Dynamics 

In this Section, we discuss the microscopic view of a dynamic behaviour of in- 
teracting many-body systems on a lattice. It was recognized for many years 
that the strong correlation in solids exist between the motions of various par- 
ticles ( electrons and ions, i.e. the fermion and boson degrees of freedom ) 
which arise from the Coulomb forces. The most interesting objects are met- 
als and their compounds. They are invariant under the translation group 
of a crystal lattice and have lattice vibrations as well as electron degrees of 
freedom. There are many evidences for the importance of many-body effects 
in these systems. Within the Landau semi-phenomenological theory it was 
suggested that the low-lying excited states of an interacting Fermi gas can be 
described in terms of a set of "independent quasi-particles". However, this 
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was a phenomenological approach and did not reveal the nature of relevant 
interactions. 



6.1 Green Function Picture of Quasi-Particles 

An alternative way of viewing quasi-particles, more general and consistent, 
is through the Green function scheme of many-body theory [Q], which we 
sketch below for completeness and for pedagogical reasons. 
We should mention that there exist a big variety of quasi-particles in many- 
body systems. At sufficiently low temperatures, few quasi-particles are ex- 
cited, and therefore this dilute quasi-particle gas is nearly a non-interacting 
gas in the sense that the quasi-particles rarely collide. The success of the 
quasi-particle concept in an interacting many-body system is particularly 
striking because of a great number of various applications. However, the 
range of validity of the quasi-particle approximation, especially for strongly 
interacting lattice systems, was not discussed properly in many cases. In 
systems like simple metals, quasi-particles constitute long-lived, weakly in- 
teracting excitations, since their intrinsic decay rate varies as the square of 
the dispersion law, thereby justifying their use as the building blocks for the 
low-lying excitation spectrum. 

Unfortunately, there are many strongly correlated systems on a lattice for 
which we do not have at present the truly the first-principles proof of a 
similar correspondence of the low-lying excited states of noninteracting and 
interacting systems, adiabatic switching on of the interaction, a simple ef- 
fective mass spectrum, long lifetimes of quasi-particles, etc. These specific 
features of strongly correlated systems are the main reason of why the usual 
perturbation theory starting from noninteracting states does not work prop- 
erly. Many other subtle nonanalytic effects which are present even in normal 
systems have the similar nature . This lack of a rigorous foundation for the 
theory of strongly interacting systems on a lattice is not only a problem 
of the mathematical perfectionism, but also that of the correct physics of 
interacting systems. 

As we mentioned earlier, to describe a quasi-particle correctly, the Green 
functions method is a very suitable and useful tool. What concerns us here 
are formal expression for the single-particle GF (38) and the correspond- 



ing quasi-particle excitation spectrum. From the equation ( 24) it is thus 
seen that the GF is completely determined by the spectral weight function 
A(lo). The spectral weight function reflects the microscopic structure of the 
system under consideration. The other term in ( is a separation of the 
purely statistical aspects of GF. From the equation ( ^0|) it follows that the 
spectral weight function can be written formally in terms of many-particle 
eigenstates. Its Fourier transform origination ( |l^ is then the density of 
states that can be reached by adding or removing a particle of a given mo- 
mentum and energy. 
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Consider a system of interacting fermions as an example. For a nonin- 
teracting system, the spectral weight function of the single-particle GF 
Gk{oj) =<< Ofco-i^feo- simple peaked structure 

. For an interacting system, the spectral function Ai.{uj) has no such a simple 
peaked structure, but it obeys the following conditions 

^fc(w) > 0; Ak{uj)duj =<[aka,alj+>=l 



Thus, we can see from these expressions that for a noninteracting system, the 
sum rule is exhausted by a single peak. A sharply peaked spectral function 
for an interacting system means a long-lived single-particle-like excitation. 
Thus, the spectral weight function was established here as the physically 
significant attribute of GF. The question of how best to extract it from a 
microscopic theory is the main aim of the present review. 
The GF for the non-interacting system is Gk{uj) = {uj — e^)^^. For a weakly 
interacting Fermi system, we have Gk{i^) = {uj — eu — Mk{uj))^'^ where Mk{uj) 
is the mass operator. Thus, for a weakly interacting system, the (^-function 
for Ak{oo) is spread into a peak of finite width due to the mass operator. We 
have 

Mk{uJ ± ie) = ReMk{uj) T ImMk{uj) = Afe(u;) ^ Tk{oj) 
The single-particle GF can be written in the form 

(76) Gk{u:) = {a; - [e^ + Afc(a;)] ± Tk{uj)}-^ 

In the weakly interacting case, we can thus find the energies of quasi-particles 
by looking for the poles of single-particle GF ([76|) 

LO = ek + Ak{uj) ± Tkiuj) 

. The dispersion relation of a quasi-particle 

e{k) = ek + Ak[e{k)]±Tk[e{k)] 

and the lifetime l/F^ then reflects the inter-particle interaction. It is easy 
to see the connection between the width of the spectral weight function and 
decay rate. We can write 

(77) Akiu;) = (exp(/5cu) + l)-\-i)[Gkiu; + ie) - Gk{u; - ie)] = 



uj-{ek + Ak{uj))]^ + Tl{uj) 



In other words, for this case, the corresponding propagator can be written 
in the form 

Gk(t) ^ eyip{-ie{k)t) exp(-rfct) 
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This form shows under which conditions, the time-development of an in- 
teracting system can be interpreted as the propagation of a quasi-particle 
with a reasonably well-defined energy and a sufficiently long lifetime. To 
demonstrate this, we consider the following conditions 

Afc[e(fc)] «e(A;); rfc[e(A:)] « e(A:) 

Then we can write 

(78) Gfc(w) = 

[^-e(fc)][l-^Uew]+^r,[e(A;)] 

where the renormalized energy of excitations is defined by 

e{k) = ek + Ak[e{k)] 

In this case, we have, instead of ( |77 



(79) Ak{uj) = 

1 1 r dAk(uj) , , 1 2T(k) 



As a result, we find 



(80) Gk it)=« aka (t) ; aL »= 

= -i9it)exp{-ie{k)t)exp{-mt)[l - 

A widely known strategy to justify this line of reasoning is the perturba- 
tion theory Q. A detailed analysis of various successful approximations for 
the determination of excited states in the framework of the quasi-particle 
concept and the Green functions method in metals, semiconductors, and 
insulators was done in review paper [ 36 1. 



There are examples of weakly interacting systems, i.g. the superconducting 
phase, which are not connected perturbatively with noninteracting systems. 
Moreover, the superconductor is a system in which the interaction between 
electrons qualitatively changes the spectrum of excitations. However, quasi- 
particles are still of use even in this case, due to the correct redefinition of 
the relevant generalized mean field which includes the anomalous averages 
(see (|72[)). In a strongly interacted system on a lattice with complex spec- 
tra, the concept of a quasi-particle needs a suitable adaptation and a careful 
examination. It is therefore useful to have the workable and efficient IGF 
method which, as we shall see, permits one to determine and correctly sep- 
arate the elastic and inelastic scattering renormalizations through a correct 
definition of the generalized mean field and to calculate real quasi-particle 
spectra, including the damping and lifetime effects. A careful analysis and 
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detailed presentations of the IGF method will provide an important step to 
the formulation of the consistent theory of strongly interacting systems and 
the justification of approximate methods presently used within equation-of- 
motion approaches. These latter remarks will not be substantiated until 
next Sections, but it is important to emphasize that the development which 
follows is not a merely formal exercise but essential for the proper and con- 
sistent theory of strongly interacting many-body systems on a lattice. 



6.2 Spin- Wave Scattering Effects in Heisenberg Ferromagnet 

In this Section, we briefly describe , mainly for pedagogical reasons, how 
the formulation of the quasi-particle picture depends in an essential way on 
an analysis of the sort introduced in Section 3.1. We consider here the most 



studied case of a Heisenberg ferromagnet 1 47] with the Hamiltonian (^) and 
the equation of motion (|66|). In an earlier discussion in Sections 4.11 and 5.3, 
we described the Tyablikov decoupling procedure (|6^) based on replacing Sf 
by < Sf > in the last term of (^) . We also discussed an alternative method 
of decoupling proposed by Callen (|69|) . Both these decoupling procedures 
retain only the elastic spin-wave scattering effects. But for our purposes, 
it is essential to retain also the inelastic scattering effects, and therefore, 
we must carefully identify and separate the elastic and inelastic spin-wave 
scattering. This is directly related with the correct definition of generalized 
mean fields. Thus, the purpose of the present consideration is to justify the 
use of IGF method for the self-consistent theory of spin-wave interactions. 
The irreducible part of GF is introduced according to the definition ( |30| ) as 

<< {S^ Sg—S'^ Si)\Sj >>=<< {SfSg—S^Si)—AigS^—AgiS'^\Sj » 
(81) 

Here the unknown quantities Aig are defined on the basis of orthogonality 
constraint (|3l| ) 



<[{stsi-s^sfY^^\s.] >=0 

We have (i / g) 

_ 2 < SfSI > + < SrS+ > 
(^2) — Agi — 2 < > 

The definition (see eq.(p3D ) of a generalized mean field GF G*^^ is given 
by the equation 



(83) ^Gf/^ = 2<S'>5^j+Y^ JigA.g {Gf/ - G- 



MF r<MF\ 

g 



From the Dyson equation in the form ( pTI) we find 

(84) Mij = {Pi,y = 

< 2S' >-2 J2 JigJij « {StS'g - S^Sff'^iiStS'g - 5+5f )(^'-))t 



gi 
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where the proper (p) part of the irreducible GF is defined by the equation 



m 



gi 

( in the diagrammatic language, this means that it has no parts connected 
by one G*^^-line) . The formal solution of the Dyson equation is of the form 
(38): 

(85) Gij{u:) = 
2<S' > A^-^ ^exp[iA;(i?i - Rj)][uj - uj{k) - 2 < S' > Mk{uj)]-^ 

k 

The spectrum of spin excitations in the generalized mean field approximation 
is given by 

(86) uik) = N-^ J2 'JigAig{l - eM^KRi " Rj)]} 

Now it is not difficult to see that the result ( p6| ) includes both the simplest 
spin- wave dispersion law ( ^T|) and the result of Tyablikov decoupling (67) 
as the limiting cases 

(87) uj{k) =< S' > (Jo - Jk) + 
(< 2S^ > iV)-i ^(J, - J,_,)(V-+ + 2Vf ) 

where 

= E < s^^t > ^MMRi - Rj)] 



It is seen that due to the correct definition of generalized mean fields (82) 
we get the spin excitation spectrum in a general way. In the hydrodynamic 
limit, it leads to Lo{k) ~ k'^. The procedure is straightforward, and the de- 
tails are left as an exercise. 

Let us remind that till now no approximation has been made. The expres- 



sions (84), (ISq), and (^) are very useful as the starting point for approx- 
imate calculation of the self-energy, a determination of which can only be 
approximate. To do this, it is first necessary to express, using the spectral 



theorem (p^), the mass operator (84) in terms of correlation functions 

< 2S' > Mkiio) = 



1 /■+°° duj' /•+00 

— / -(exp(/3a;') - 1) / dtexp{iLo't) 

27r J-oo 00 — UJ J-oo 



1— < ((5+(t)5|(t) - 5+(t)5f (t))("'))t|(S+5,^ - )(-) >(P) 



< 2S^ > 
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This representation is exact, and only the algebraic properties were used to 
derive it. Thus, the expression for the analytic structure of the single-particle 
GF ( or the propagator ) can be deduced without any approximation. A 
characteristic feature of eq.(|8^) is that it involves the higher-order GFs. 
A whole hierarchy of equations involving higher-order GFs could thus be 
rewritten compactly. Moreover, it not only gives a convenient alternative 
representation, but avoids some of the algebraic complexities of higher-order 
Green-function theories. Objective of the present consideration is to give a 
plausible self-consistent scheme of the approximate calculation of the self- 
energy within the IGF method. To this end, we should express the higher- 
order GFs in terms of the initial ones, i.e. find the relevant approximate 
functional form 

M w F[G] 

It is clear that this can be done in many ways. As a start, let us consider how 
to express higher-order correlation function in (^8|) in terms of the low-order 
ones. We use the following form! 



(89) < - 5;(t)5f (t))(-))t|(5+5| - >(P)^ 

We find 

(90) < 2S' > Mk{Lo) = 

1 /■+°° duj' f~^^ 

-{exp{f3io') - 1) dteyip{iuj't) 



27r j_oo uj — uj' 



-±—(^^^.^^(t)^-+(t) - vf/(t)^-+(t) - r,tmr,^{t) + i^frm-;{t)) 

It is reasonable to approximate the longitudinal correlation function by its 
static value ipjfit) ~ (0). The transversal spin correlation functions are 
given by the expression 

(91) i^Ji^it) = 

/OO 
— [exp(/5u;) - 1]"^ exp(iu;t)(-2/m << S:^\Sj' »u,+ie) 



After the substitution of eq.( |9l|) into eq.( for the self-energy, we find 
an approximate expression in the self-consistent form, which, together with 
the exact Dyson equation (^5|), constitute a self-consistent system of equa- 
tions for the calculation of the GF. As an example, we start the calculation 
procedure ( which can be made iterative ) with the simplest first "trial" 
expression 

{-2Im « S+\Sr »^+,,) « 6{uj - oj{k)) 
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After some algebraic transformations we find 

(92) < 2S' > Mk{u) ^ iV-i ^(J, - Jk-q?{oJ - uj{q - fe))-^^ 

This expression gives a compact representation for the self-energy of the 
spin-wave propagator in a Heisenberg ferromagnet. The above calculations 
show that the inelastic spin-wave scattering effects influence the single- 
particle spin-wave excitation energy 

oj{k, T) = Lo{k) + ReMk{uj{k)) 

and the energy width 

Tk{T)=ImMk{io{k)) 

Both these quantities are observable, in principle, via the ferromagnetic reso- 
nance or inelastic scattering of neutrons. There is no time to go into details 
of this aspect of spin-wave interaction effects. It is worthy to note only 
that it is well known that spin-wave interactions in ferromagnetic insulators 
have a relatively well-established theoretical foundation, in contrast to the 
situation with antiferromagnets. 



7 Heisenberg Antiferromagnet at Finite Temper- 
atures 

As it is mentioned above, in this article, we describe the foundation of 
the IGF method, which is based on the equation-of-motion approach. The 
strength of this approach lies in its flexibility and applicability to systems 
with complex spectra and strong interaction. The microscopic theory of 
the Heisenberg antiferromagnet is of great interest from the point of view 
of application to any novel many-body technique. This is not only because 
of the interesting nature of the phenomenon itself but also because of the 
intrinsic difficulty of solving the problem self-consistently in a wide range of 
temperatures. In this Section, we briefly describe how the generalized mean 
fields should be constructed for the case of the Heisenberg antiferromagnet, 
which become very complicated when one uses other many-body methods, 



like the diagrammatic technique |87]. Within our IGF scheme, however, the 



calculations of quasi-particle spectra seem feasible and very compact. 



7.1 Hamiltonian of the Model 

The problem to be considered is the many-body quasi-particle dynamics of 
the system described by the Hamiltonian [46| 

(93) h = -Wy1 "^""'(^ - 3%J,a' = - J E E Jr's,J-,c.' 



ij aa' 



1 cxa' 
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This is the Heisenberg-Neel model of an isotropic two-sublattice antiferro- 
magnet (the notation is shghtly more general than in Section 4.1.2 ). Here 
Sia is a spin operator situated on site i of sublattice a, and J"" {i — j) is the 
exchange energy between atoms on sites Ria and Rja'] ol^ a' takes two values 
(a, b) . It is assumed that all of the atoms on sublattice a are identical, with 
spin magnitude Sa- It should be noted that, in principle, no restrictions are 
placed in the Hamiltonian ( p3[ ) on the number of sublattices, or the number 
of sites on a sublattice. What is important is that sublattices are to be dis- 
tinguished on the basis of differences in local magnetic characteristics rather 
than merely differences in geometrical or chemical characteristics. 
Let us introduce the spin operators = Sf^±iSf^. Then the commutation 
rules for spin operators are 

For an antiferromagnet, an exact ground state is not known. Neel [^] in- 
troduced the model concept of two mutually interpenetrating sublattices 
to explain the behaviour of the susceptibility of antiferromagnets. How- 
ever, the ground state in the form of two sublattices ( the Neel state ) is 
only a classical approximation. In contrast to ferromagnets, in which the 
mean molecular field is approximated relatively reasonably by a function 
homogeneous and proportional to the magnetisation, in ferri- and antifer- 
romagnets, the mean molecular field is strongly inhomogeneous. The local 



molecular field of Neel |67] is a more general concept. Here, we present the 
calculations of the quasi-particle spectrum and damping of a Heisenberg 
antiferromagnet in the framework of the IGF method. 
In what follows, it is convenient to rewrite ( |93| ) in the form 

(94) H=-WY: ir'iS^aSZ,^' + 



where 



= i/2( + r 



It will be shown that the use of "anomalous averages" which fix the Neel 
vacuum makes it possible to determine uniquely generalized mean fields 
and to calculate, in a very compact manner, the spectrum of spin-wave 
excitations and their damping due to inelastic magnon-magnon scattering 
processes. A transformation from the spin operators to Bose (or Pauli ) 
operators is not required. 



7.2 Quasi-Particle Dynamics of Heisenberg Antiferromagnet 

In this section, to make the discussion more concrete, we consider the re- 
tarded GF of localized spins defined as G^^{t - t') =« A{t),B{t') » . 
Our attention is focused on the spin dynamics of the model. To describe 
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the spin dynamics of the model ( self-consistently, one should take into 
account the full algebra of relevant operators of the suitable "spin modes" 
( "relevant degrees of freedom" ) which are appropriate for the case. This 

relevant algebra should be described by the 'spinor' A = { '+), B = A^ , 

^kb 

according to the IGF strategy of Section 3. 

Once this has been done, we must introduce the generalized matrix GF of 
the form 

(95) (^^ S'fl'' = Gik;u;) 

\«Stb\S^ka» «Stb\S-kb»J 

To show the advantages of the IGF in the most full form, we carry out the 
calculations in the matrix form. 

To demonstrate the utility of the IGF method, we consider the following 
steps in a more detailed form. Differentiating the GF << S^JB » with 
respect to the first time, t, we find 



(96) '^«^.+J 




2<5,^>1 1 



yl/o ^kql^ab >>uj 

g 

'^lkj2 Yl C « Sll\Bab »u 



where — (S^.g ^Sgj, 5'g^5'^_g,a)- 
In (96), we introduced the notation 



Let us define the irreducible (ir) operators as (equivalently, it is possible to 
define the irreducible GFs) 

(97) (Sf^) = Sf^ - AfSt, + A\%St, 

(98) {Sl^) ^^^^ = Sl^-N^I^<S^^>5,^o 

The choice of the irreducible parts is uniquely determined by the "orthogo- 
nality" constraint (31) 

(99) < [{St\) I '^-^'^ 




From eq.(|99|) we find that 

2< (5i,j("')(S^\)(-) > + < > 



(100) A; 



« 2iVV2 < Sf, > 
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By using the definition of tlie irreducible parts ( |97| ) , the equation of motion 
(96) can be exactly transformed to the following form 

(101) (W - OJaa) « S^jBab »u> +^ab « S^^\Bab »u> = 

(102) {UJ - UJbb) « S^^\Bba »u> +OJba « S^jBba »u> = 

The following notation was used: 

(103) u:aa = ((C - ID < SI > <SI> + 



(104) u;ab={lf<S^,>+J2lt^A 



ba 
Nq 

q 



(105) Af^ = N'^/^Af 

(106) $fnfe) = 

1 'y=a,b 

To calculate the irreducible GFs on the right-hand sides of eqs. ( |101| ) and 
(|102| ), we use the device of differentiating with respect to the second time t'. 
After introduction of the corresponding irreducible parts into the resulting 
equations, the system of equations can be represented in the matrix form 
which can be identically transformed to the standard form (jsj) 

(107) G{k, Lo) = Go{k, uj) + Go{k, u)P{k, uj)Go{k, u) 

Here we introduced the generalized mean- field (GMF) GF Go and the scat- 
tering operator P according to the following definitions 

(108) Go = n^^i 



(109) P = 
1 (« )(fc)|$r)^(fc) » « $f )(A;)|<I>f )t(A;) »' 

4 < 5,- >2 c^f )(A;)|$^)t(A;) » « $f )(A;)|<I>f )^(A;) », 

where 

(110) C^=((^-^aa) U^ab 

\ ^ab {l^ - i^bb) 
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The Dyson equation can be written exactly in the form (^) where the mass 
operator M is of the form 

(111) M{k,oj) = [P{k,uj))^P'^ 

It follows from the Dyson equation that 

P{k,uj) = M{k,uj) + M{k,u;)Go{k,Lu)P{k,u;) 

Thus, on the basis of these relations, we can speak of the mass operator M 
as the proper part of the operator P by analogy with the diagram technique, 
in which the mass operator is the connected part of the scattering operator. 
As it is shown in Section 3, the formal solution of the Dyson equation is of 
the form (38). Hence, the determination of the full GF G was reduced to 
the determination of Gq and M. 



7.3 Generalized Mean-Field GF 



From the definition (108), the GF matrix in the generalized mean- field ap- 
proximation reads 



(112) 



Gr{k,io) Gf{k,u;) 
G(j"(A;,w) G','ik,u;) 



Go-- 

< Sa> / {UJ - Uaa) ^ab 

detn \ ^ab {(^ - ^^bb) 



where 



detO. = {u) - U}aa){i^ - (^bb) - i^aa^ab 

We find the poles of GF ( |112| ) from the equation 

detVt = 

from which it follows that 

(113) ^±(A:) = ±^{u:lM-^lbm 

It is convenient to adopt here the Bogoliubov (n, t>)-transformation notation 
by analogy with that of Section 4.1.2. The elements of the matrix GF 
Go(A;,a;) are found to be 

u'^{k) v'^ik) 



(114) Gr{k,u;)=2<S^,> 



(115) G^'{k,u;) 
where 



2 < 5,^ > 



cj— a;-l_(A;) uj — uj^{k)- 
-u{k)v{k) ^ u{k)v{k) 



uj — uj^{k) uj — uj-{k) 



--G',-{k,u;) 



u\k) = 1/2[(1 - 7fc')-V2 + 1]. ^2(^) ^ i/2[(l - ^l)-'/' - 1] 
(116) 7fc = -Eexp(iA^i?^); /,"'^ = lf = 
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The simplest assumption is that each sublattice is s.c. and ujaaik) = 
{a = a,b). Although that we work in the GFs formalism, our ex- 
pressions ( |114[ ), (115) are in accordance with the results of the Bogoliubov 



(u,v)-transformation , but, of course, the present derivation is more general. 
However, it is possible to say that we diagonalized the generalized mean-field 
GF by introducing a new set of operators. We used the notation 



(117) 



S+{k) = UkS+^ + ?;fc5+ ; S+{k) = + ^^^5+ 



This notation permits us to write down the results in a compact and conve- 
nient form, but all calculations can be done in the initial notation too. 
The spectrum of elementary excitations in the GMF approximation for an 
arbitrary spin S is of the form 

1 ■ afe 



(118) ^(fc) = iz<si>[i- ^ ^ 7.<J - rt) 

where Iq = zl'jq, and z is the number of nearest neighbors in the lattice. The 
first term in (118) corresponds to the Tyablikov approximation ( c/.(H^). 



The second term in (|11^ ) describes the elastic scattering of the spin-wave 
quasi-particles. At low temperatures, the fluctuations of the longitudinal 
spin components are small, and, therefore, for ( |118| ) we obtain 



(119) u;{k) ^ ISz[l - C{T)]^{l--fl) 

The function C(T) determines the temperature dependence of the spin- wave 
spectrum 

(120) C{T) = ^ Y.(< S-,aSta > +7, < SZ,,5+ >) 

In the case when C(T) — > 0, we obtain the result of the Tyablikov decoupling 
for the spectrum of the antiferromagnons 

(121) cvik) ^I<Sl> z^{l-^l) 

In the hydrodynamic limit, when uj{k) ~ D{T)\k\, we can conclude that 
the stiffness constant D{T) = zIS{\ — C(T)) for an antiferromagnet de- 
creases with temperature because of the elastic magnon-magnon scattering 
as T^. To estimate the contribution of the inelastic scattering processes, it 
is necessary to take into account the corrections due to the mass operator. 

7.4 Damping of Quasi-Particle Excitations 

An antiferromagnet is a system with a complicated quasi-particle spectrum. 
The calculation of the damping due to inelastic scattering processes in a 
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system of that sort has some important aspects. When calculating the 
damping, it is necessary to take into account the contributions from all 
matrix elements of the mass operator M 

M = Gq^ - G-^ 

It is then convenient to use the representation in which the generalized mean 
field GF has a diagonal form . In terms of the new operators 5i and 52, the 
GF G takes the form 

r(k.>,\-(«~^tm~Sn-k)» «5+(fe)|52-(-fc)»\ _ /Gn Gi2 
u{K,uj) [^^^s+ik)\S^i-k)» «S+ik)\S^i-k)»J \G2i G22 

In other words, the damping of the quasi-particle excitations is determined 
on the basis of a GF of the form 

(122) Gn{k,io)- < > 



oj-Lo{h) - 2<S^> T.{k,uj) 
Here, the self-energy operator T,{k,uj) is determined by the expression 

(123, E(...) = M„(...)- ^;|r^<S^ 

In the case when A;, u; — > 0, one can be restricted to the approximation 

(124) S(A:,a;) « Mu{k,io) = ulMaa + VkUk{Mab + Mba) + vlUbb 

It follows from ( |111D that to calculate the damping, it is necessary to find the 
GFs « ^^j\k)\^^'^Hk) 

>>. As an example, we consider the calculation 
of one of them. By means of the spectral theorem (pTj), we can express the 
GF in terms of the correlation function < ^a'^^\k)^^^\k,t) >. We have 

(125) « «>f-)(A:)|«>i*''-)t(A;) »= 

1 z'+cxD fj, J r-\-oo 

— / ;(exp(/3a;') - 1) / dtexp(ia;'t) < <I>f )t(A;)cI>f )(fc, t) > 

27r J -00 uj — uj' J -00 

Thus, it is necessary to find a workable "trial" approximation for the corre- 
lation function on the r.h.s. of ( |125| ). We consider an approximation of the 
following form 



(126) <(5i,.)(^'-)5:(._,0a5;-,0a(^)(5.%W)^^^^ 



4iV52 p 



where ipqtbi^) ~< ^-qa^qbi^) >• analogy with the diagram technique, 
we can say that the approximation (|126|) corresponds to the neglect of the 
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vertex corrections to the magnon-magnon inelastic collisions. Using ( |126| ) 
in ( |125| ), we obtain 

(127) « $i^'-)(A:)|«>i*'-)t(fc) »~ 

1 ^ f dujiduj2duJs 

16NS^ ^ J ll> - uJi - UJ2 + 0J3 

[--ImGaaik - q,uJi)][--ImGbbiq + P,uJ2)][--ImGbb{p,uJ3)] 

TT TT TT 

where 

(128) F{uJuuj2,uJ3) = N{uj2)[N{lo3) - N{uJi)] + [1 + N{uji)]N{lo3) 



Equations (p7\), (111), and ( p.27D constitute a self-consistent system of equa- 
tions. To solve this system of equations, we can, in principle, use any con- 
venient initial representation for the GF, substituting it into the right-hand 



side of eq. (127). The system can then be solved iteratively. To estimate 
the damping, it is usually sufficient, as the first iteration, to use the simplest 
single-pole approximation 

(129) - -ImG{k, uj) « 6{uj - uj{k)) 

TT 

As a result, for the damping of the spin- wave excitations we obtain 
(130) r(A;, u) = -2SImT.{k, u) = 

Np{l + iVg+p)(l + Nk-g)Muik,p; k-q,p + q)5{uj - uj{k - q) + i^{p)) 

qp 

The explicit expression for Mn is given in ref. [^ ]. In our approach, it 
is possible to take into account the inelastic scattering of spin waves due 
to scattering by the longitudinal spin fluctuations too |3^. In general, the 
correct estimates of the temperature dependence of the damping of antifer- 
romagnons depend strongly on the reduced temperature and energy scales 
and are rather a nontrivial task. However, under the normal conditions, 
the damping is weak u}{k)/T ~ 10^ — 10'^, and the antiferromagnons are the 



well-defined quasi-particle excitations 1 88]. 



In summary, in this Section, we have shown that the IGF method permits 
us to calculate the spectrum and the damping for a two-sublattice Heisen- 
berg antiferromagnet in a wide range of temperatures in a compact and 
self-consistent way. At the same time, a certain advantage is that all the 
calculation can be made in the representation of spin operators for an arbi- 
trary spin S. The theory we have developed can be directly extended to the 
case of a large number of magnetic sublattices with inequivalent spins, i.e.. 
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it can be used to describe the complex ferrimagnets. 

In the framework of our IGF approach, it was shown that the mean fields 
in an antiferromagnet must include the "anomalous" averages which rep- 
resent the local nature of the Neel molecular fields. Thus, the mean field 
in an antiferromagnet, like the mean field in a superconductor, has a more 
complicated structure. 



8 Quasi-Particle Dynamics of Lattice Fermion Mod- 
els 

8.1 Hubbard Model. Weak Correlation 

The concept of GMFs and the relevant algebra of operators from which GFs 
are constructed are important for our treatment of electron correlations in 
solids. It is convenient (and much shorter) to discuss these concepts for 
weakly and strongly correlated cases separately. First, we should construct 



a suitable state vector space of a many-body system [|10[ . The fundamental 
assumption implies that states of a system of interacting particles can be 
expanded in terms of states of non-interacting particles This approach 
originates from perturbation theory and finds support for weakly interacting 
many-particle systems. For the strongly correlated case, this approach needs 
a suitable reformulation, and just at this point, the right definition of the 
GMFs is vital. Let us consider the weakly correlated Hubbard model (^9[). In 
some respect, this case is similar to the ordinary interacting electron gas but 
with very local singular interaction. The difference is in the lattice ( Wannier 
) character of electron states. It is shown below that the usual creation 
and annihilation ajo- second-quantized operators with the properties 

^1^(0)^^(1). = M/(0) 

ai^'(0)=0; aj^S'^ = (i ^ j) 

are suitable variables for description of a system under consideration. Here 
^'(0) and ^-(1) are vacuum and single-particle states, respectively. The ques- 
tion now is how to describe our system in terms of quasi-particles. For a 
translationally invariant system, to describe the low-lying excitations of a 
system in terms of quasi-particles [Q] , one has to choose eigenstates such that 
they all correspond to a definite momentum. For the single-band Hubbard 
model (|49|), the exact transformation reads 

aj:^ = N'^/"^ ^ exp(-iA;^i)ai^ 

i 

Note that for a degenerate band model, a more general transformation is 
necessary 1 91]. Then the Hubbard Hamiltonian (^) in the Bloch vector state 
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space is given by 

(131) H = J2 e{k)alak^ + U/2N ^ 

ka pqrs 

If the interaction is weak, the algebra of relevant operators is very simple: it 
is an algebra of a non-interacting fermion system (flfco-i alo-' ""fco- ~ (^ka^ka)- 
To calculate of the electron quasi-particle spectrum of the Hubbard model 
in this limit, let us consider the single-electron GF defined as 

(132) Gk^t - t') =« aka, aL »= "^^(^ " < Mt), aL(*')]+ > 

The equation of motion for the Fourier transform of GF Gkai^) is of the 
form 



(133) {u;-ek)GkA^) = l + U/NY^ 



pq 

Let US introduce an "irreducible" GF in the following way 

(ir) t it _ 

0'k+paO"p-\-q—a-^(]~o-\^k(7 ^^(^ — 

(134) << afc+pcrap_,_g_^ag_o-|aio- >><^ ~^p,o < ^q-a > Gka 



The irreducible (ir) GF in ( |134D is defined so that it cannot be reduced to 
GF of lower order with respect to the number of fermion operators by an 
arbitrary pairing of operators or, in other words, by any kind of decoupling. 
Substituting ( [L34D into (|133| ), we obtain 



(135) GU-^) = Gt^iu;) + 

pq 

Here we introduced the notation 

(136) Gl'fiuj) = {iv- eika))-'; e{ka) = e{k) + U/NJ2< Uq., > 

In this paper, for brevity, we confine ourselves to considering the paramag- 
netic solutions, i.e. < n„ >=< >. To calculate the higher-order GF 
on the r.h.s. of ( |135| ), we have to write the equation of motion obtained by 
means of differentiation with respect to the second variable t'. Constraint 
(31) allows us to remove the inhomogeneous term from this equation for 

^''U<Ait),al{t')». 

For the Fourier components, we have 

{u - 6(A;))(-) « A\al ».=<(^^) [^,4j+ > + 

(137) U/NJ2^'^^ « A\al_^ar+s-aal^,^ »^ . 

rs 
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The anticommutator in ( |137| ) is calculated on the basis of the definition of 
the irreducible part 

(138) < [ak+paO-p+q-a^q-rr— < 0.p+q~a^q~<T > '^fc+pcri ^ 

If one introduces the irreducible part for the r.h.s. operators by analogy with 
expression ( |134D , the equation of motion (133) takes the following exact form 
( cf. eq.®) 



(139) GU^) = Gt^{uj) + Gf/{u:)PUu)Gt'f{uj) 

where we introduced the following notation for the operator P (35) 



2 



(140) Pfc.(^) = ^^Z)£)(p,(?|r,.,; 



Ar2 

pqrs 



({ir) t it t 

2^ { '^k+pa'J'p+q-a^q-t^l^r-a^r+s-aO'k+sa >>a; J 



pqrs 



To define the self-energy operator according to ( |3q ) , one should separate the 
"proper" part in the following way 



Dia\p^l\r,s;uj) = l''^^\p, q\r, s;u;) 
(141) +^ E Lt\p,q\r's'-,u.)GZ^{u;)D^-J{p',q'\r,s-,u) 

r's'p'q' 



Here L^i^^\p,q\r,s;io) is the "proper" part of GF D^i!^\p,q\r,s;io) which, in 
accordance with the definition (36), cannot be reduced to the lower-order 
one by any type of decoupling. We find 

(142) Gka = Gffiio) + Gff{Lo)Mk.{Lo)GkA^) 

Equation ( |142| ) is the Dyson equation for the single-particle double-time 
thermal GF. According to (38), it has the formal solution 

(143) Gka{oo) = [oJ- e{ka) - Mka{uj)]-^ 
where the self-energy operator M is given by 

(144) MU^) = ^ E LfJ{p,q\r,s;u:) = 



Ar2 

pqrs 



C/2 
2 

pqrs 



We wrote explicitly equation (|14C| ) for P and equation (144) for M to illus- 
trate the general arguments of Section 3 and to give concrete equations for 
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determining both the quantities, P and M. 

The latter expression ( |144| ) is an exact representation ( no decoupling was 
made till now ) for the self-energy in terms of higher-order GF up to sec- 
ond order in U { for the consideration of higher-order equations of motion, 
see ref. p5| ). The explicit difference between P and M follows from the 
functional form (|38|). Thus, in contrast to the standard equation-of- motion 
approach, the calculation of full GF was substituted by the calculation of 
the mean-field GF G^^^ and the self-energy operator M. The main rea- 
son for this method of calculation is that the decoupling is only introduced 
into the self-energy operator, as it will be shown in a detailed form below. 



The formal solution of the Dyson equation (38) determines the right refer- 
ence frame for the formation of the quasi-particle spectrum due to its own 
correct functional structure. In the standard equation-of-motion approach, 
that structure could be lost by using decoupling approximations before ar- 
riving at the correct functional structure of the formal solution of the Dyson 
equation. This is a crucial point of the IGF method. 

The energies of electron states in the mean-field approximation are given 
by the poles of . Now let us consider the damping effects and finite 

lifetimes. To find an explicit expression for the self-energy M ( |144D , we have 
to evaluate approximately the higher-order GF in it. It will be shown below 
that the IGF method permits one to derive the damping in a self-consistent 
way simply and much more generally than within other formulations. First, 
it is convenient to write down the GF in ( |144 ) in terms of correlation func- 
tions by using the spectral theorem (26) 

— / -{expif3ij') + l) expiiu't) 

(145) < al^^^{t)al._^{t)ar+s~a{t)ak+paal+g-aaq-a > 

Further insight is gained if we select the suitable relevant "trial" approxi- 
mation for the correlation function on the r.h.s. of ( |145| ). In this paper, we 
show that the earlier formulations based on the decoupling or/and diagram- 
matic methods can be obtained from our technique but in a self-consistent 
way. It is clear that a relevant trial approximation for the correlation func- 
tion in ( |145D can be chosen in many ways. For example, the reasonable and 
workable one can be the following "pair approximation" that is especially 
suitable for a low density of quasi-particles: 

(146) 

^k+s,k+p^r,q^r+s,p+q 

Using (|146|) and (|145D in (|144|) we obtain the self-consistent approximate 



expression for the self-energy operator ( the self-consistency means that we 



53 



express approximately the self-energy operator in terms of the initial GF, 
and, in principle, one can obtain the required solution by a suitable iteration 
procedure ): 



(147) 



El 

Ar2 



E 

pq 



duJiduj2duJ3 

U> + UJi — UJ2 — UJs 



n{uJ2)n{iO^) + n(a;i)(l - n{uJ2) - ^(tJs)) 9p+q-a{'~^l)9k+pa{l^2)gq^a{l^-i) 



where we used the notation 
1 



TT 



ImGka{'^ + i^)] n{uj) = [exp(/?L<;) + 1] 



Equations ( |147D and ( |142| ) constitute a closed self-consistent system of equa- 
tions for the single-electron GF of the Hubbard model in the weakly corre- 



lated limit. In principle, we can use, on the r.h.s. of (147), any workable first 
iteration-step form of the GF and find a solution by iteration ( see Appendix 
D ). It is most convenient to choose, as the first iteration step, the following 
simple one-pole approximation: 



(148) 



Qkai^) ~ - e(fco-)) 



Then, using (p8|) in ( PItI) , we get, for the self-energy, the explicit and 
compact expression 



(149) MU^) = 1-^Y.^ 



p+q- 



-cr(l T^k+pa '^q—a) ~\~ '^k+pa'^q—a 



Ar2 



pq 



uj + e{p + qa) - e{k + pa) - e{qa) 



Formula ( |149D for the self-energy operator shows the role of correlation 
effects ( inelastic scattering processes ) in the formation of quasi-particle 
spectrum of the Hubbard model. This formula can be derived by several 
different methods, including perturbation theory. Here we derived it from 
our IGF formalism as a known limiting case. The numerical calculations of 
the typical behaviour of real and imaginary parts of the self-energy ( |149| ) 
were performed |91],|31] for the model density of states of the FCC lattice. 
These calculations and many other (see e.(7.[^,[92|,|Q) show clearly that 



the conventional one-electron approximation of the band theory is not al- 
ways a sufficiently good approximation for transition metals like nickel. A 
more concrete discussion of the numerical calculations and their comparison 
with experiments deserve a separate consideration and will be considered 
elsewhere (for a detailed recent discussion, see p6| ). 

Although the solution deduced above is a good evidence for the efficiency of 
the IGF formalism, there is one more stringent test of the method that we 



54 



can perform. It is instructive to examine other types of possible trial solu- 
tions for the six-operator correlation function in the eq. (|145D . The approx- 
imation we propose now reflects the interference between the one-particle 
branch of the spectrum and the collective ones: 

< al_^^^^{t)al_„{t)ar+s-a{t)ak+paal+g^„aq^a >^'''^f« 

(150) < al._^{t)aq^^ >< al^^^{t)ar+s-a{t)ak+paal+q-a > 

It is seen that the three contributions in this trial solution describe the 
self-energy corrections that take into account the collective motions of elec- 
tron density, the spin density and the density of "doubles", respectively. 
An essential feature of this approximation is that a correct calculation of 
the single-electron quasi-particle spectra with damping requires a suitable 
incorporation of the influence of collective degrees of freedom on the single- 
particle ones. The most interesting contribution comes from spin degrees of 
freedom, since the correlated systems are often magnetic or have very well 
developed magnetic fluctuations. 

We follow the above steps and calculate the self-energy operator ( |144| ) as 

Mka{(^) = -ITT- / duJlduJ2 

iV J-oo LU - UJl — UJ2 

VexpHfc(i?i - R^)]{--Im « Sf\Sf 

(151) {--Im « ai_<^|at >>^J 

vr 

where the following notation was used: 

S+ = al^aii; S- = aj^ai^ 
It is possible to rewrite ( |151| ) in a more convenient way 

(152) MU^) = 

/■ , LO-UJ' LU\ 

—— > dio cot — h tan — - 

[-^Imx^^{k -q,uj - uj')gga{uj')^ 

Equations ( |152| ) and ( |142| ) constitute again another self-consistent system 
of equations for the single-particle GF of the Hubbard model. Note that 
both the expressions for the self-energy depend on the quasi-momentum; 
in other words, the approximate procedure does not break the momentum 
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conservation law. The fundamental importance of equations (152) and (147) 
can be appreciated by examining the problem of the definition of the Fermi 
surface. It is rather clear, because the poles uj{k,a) = e{k,a) — iTf^ of GF 
(|143| ) are determined by the equation 

u - e{ka) - Re[Mka{oj)] = 

It can be shown quite generally that the Luttinger's definition of the true 
Fermi surface Q is valid in the framework of the present theory. It is 
worthy to note that for electrons in a crystal where there is a band index, 
and a quasi-momentum, the definition of the Fermi surface is a little more 
complicated than the single-band one. Before the single particle energies 
and Fermi surface are known, one should carry out a diagonalization in the 
band index. 



8.2 Hubbard Model. Strong Correlation 

Being convinced that the IGF method can be applied successfully to the 
weakly correlated Hubbard model, we now show that the IFG approach can 
be extended to the case of an arbitrarily strong but finite interaction. This 
development incorporates main advantages of the IGF scheme and proves 
its efficiency and flexibility. 

When studying the electron quasi-particle spectrum of strongly correlated 
systems, one should take care of at least three facts of major importance: 

(i) The ground state is reconstructed radically as compared with the 
weakly correlated case. This fact makes it necessary to redefine single- 
particle states. Due to the strong correlation, the initial algebra of op- 
erators is transformed into the new algebra of complicated operators. 
In principle, in terms of the new operators, the initial Hamiltonian can 
be rewritten as a bilinear form, and the generalized Wick theorem can 
be formulated. It is very important to stress that the transformation 
to the new algebra of relevant operators reflects some important inter- 
nal symmetries of the problem, and nowadays, this way of thinking is 
formulating in the elegant and very powerful technique of the classifi- 
cation of the integrable models and exactly soluble models (c/.p^). 

(ii) The single-electron GF that describes dynamic properties, should have 
the two-pole functional structure, which gives in the atomic limit, when 
the hopping integral tends to zero, the exact two-level atomic solution. 

(iii) The GMFs have, in the general case, a very non-trivial structure. 
The GMFs functional, as a rule, cannot be expressed in terms of the 
functional of the mean particle densities. 



56 



In this section, we consider the case of a large but finite Coulomb repulsion 
U in the Hubbard Hamiltonian (|49| ) . Let us consider the single-particle GF 
( |132| ) in the Wannier basis 

(153) G^j^it - t') =« ai^it); a],(t') » 

It is convenient to introduce the new set of relevant operators []55|| 

diaa = n°'_^aia, {a = ±); n+ = riia, n~ = (1 - nj^); 

(154) ^ <^ = 1; <^nf^ = 5„/3<^; ^ di^„ = a,^ 



The new operators diaa and dj^^ have complicated commutation rules, 
namely. 

The convenience of the new operators follows immediately if one writes down 
the equation of motion for them 



(155) 



It is possible to interpret [^], |5^ both contributions to this equation as 
alloy analogy and resonance broadening corrections. Using the new operator 
algebra, it is possible identically rewrite GF ( |153| ) in the following way 

(156) Gi.aiu;) = E « d^a.\d]f,^ »^= J2 F^f^iu;) 

The equation of motion for the auxiliary matrix GF 

« di+a\d]j^^»ui « di+a\d\„„»uj\ 



(157) Ft^i[u) 



«di-a\d]^^»^ «d 



is of the following form 

(158) (EFijaiuj) - I6ij)af3 = << n'i-aO'l<T + aaia^j«-a l4/3<7 >>^ 

where the following matrix notations was used: 



(159) 



E 



{u: - E+) 








n 



— cr 









n_ 



In accordance with the general method of Section 3, we introduce by defi- 
nition the matrix IGF: 



<< Z2l\d]^^ >>, 



<< Zi2\d]_^ » 

« Z22\d\^„ » 



(160) 



A -a' 



L ija ijcr J 
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Here the notation was used: 

■Z'li = Z\2 = n'l^^ai^ + airjbii^^; Z21 = Z22 = n^_^ai„ — aiahu^^ 

It is to be emphasized that the definition ( |158| ) is the most important and 
crucial point of the whole our approach to description of the strong cor- 
relation. The coefficients A and B are determined by the orthogonality 
constraint (31), namely, 

(161) <[(Dg),^,4^J+>=0 
After some algebra, we obtain from ( |161| ) {i 7^ j) 

= a(< d\p_^ai^„ > + < di_/3_^a|"_^ >)(n^^)"^ 
[Bii]a(} = [< nf_^n"_^ > +a/3(< aiaa\_^ai-cra\„ > - 

(162) < a^„ai-,a\^yi^ >)](n^J"' 

As previously, we introduce now GMF GF F^^^; however, as it is clear from 
(|162| ), the actual definition of the GMF GF is very nontrivial. After the 
Fourier transformation, we get 

Fka^ Fk^J ah-cd\nt,c nZ^a , 
The coefficients a, 6, c, d are equal to 



(163) > 



a 
h 



(u;-E±- J2 <P)[A^^{-P) - B^^ip - q)]) 



(164) ^ = Ar-i^e(p)[^T±(_p) _ B^^{p-q)] 



Then, using the definition ( 157D , we find the final expression for GMF GF 



,MFr.\_ u - {nZ^E^ + n_^E+) - X{k) 



{UJ-E+- nZMk)){uJ -E_- ntM^)) - nZ,ntMk)^i{k) 
(165) 

Here we introduced the following notation: 
(166) ^■'^,=^Y.<P)[^^^i-P)-B^^iP-k)^ 



cr p 



(167) ^^Y.<PM^^i-P)-B^^iP-k)] 



-cr p 

X{k) = {nZ^fiXi + A3) + (nlj2(A2 + A4) 
From the equation (|165D it is obvious that our two-pole solution is more 



general than the "Hubbard HI" [^] solution and the RothpOl solution. Our 
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solution has the correct nonlocal structure and, thus, takes into account 
the non-diagonal scattering matrix elements more accurately. Those ma- 
trix elements describe the virtual "recombination" processes and reflect the 
extremely complicated structure of single-particle states which virtually in- 
clude a great number of intermediate scattering processes. 
The spectrum of mean-field quasi-particle excitations follows from the poles 
of the GF ( [165D and consists of two branches 



(168) J^{k) 



1/2[{E+ -E^+ai + hi) ± ^{E+ + E^-ai- 6i)2 - 4cd] 

where ai = uj — E± — a; bi = uj — E± — b. Thus, the spectral weight function 
Ai:^{u) of GF ( [165D consists of two peaks separated by the distance 

(169) iOi-co2 = \/{U -ai- 5i)2 - cd ^(1 - ~ ) + 0(7) 



For a deeper insight into the functional structure of the solution (165) and 



to compare with other solutions, we rewrite ( |165| ) in the following form 

db~^c\ — l d ( b da~^c\ — l 




from which we obtain for G^^^{k,uj 

MFi , _ nl^(l + c6-i) , nZA^ + da- ^' 
^^'^ ^""^ - a-db-^c + b-ca-^d ~ 

(171) ^ + ^ 

where 

(172) nt^nZ^Wt^ = N'^ ^ exp[-iA;(i2, - R^)] 

((< > + < ai-anf^a]_^ >) + 

(< nf_^nf_^ > + < ai^al^^aj-^a]^ > - < ai„ai-„a\_^a\^ >] 

are the shifts for upper and lower splitted subbands due to the elastic scat- 
tering of carriers in the Generalized Mean Field. The quantities are 
functionals of the GMF. The most important feature of the present solution 
of the strongly correlated Hubbard model is a very nontrivial structure of 



the mean-field renormalizations (171), which is crucial for understanding the 



physics of strongly correlated systems. It is important to emphasize that 
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just this complicated form of GMF is only relevant to the essence of the 
physics under consideration. The attempts to reduce the functional of GMF 
to a simpler functional of the average density of electrons are incorrect from 
the point of view of real physics of strongly correlated systems. This physics 
clearly shows that the mean-field renormalizations cannot be expressed as 
functionals of the electron mean density. To explain this statement, let us 
derive the "Hubbard I" solution Q (|5|) from our GMF solution ([l65| ). If 



we approximate ( |171| ) as 

(173) nt^nZ^W^ik) « N-'Y.i,,exp[-ikiR, - R,)] < nf^^nt^ > 
and make the additional approximation, namely. 



then solution ( 165 ) turns into the "Hubbard I" solution (|54|). This solution, 
as it is well known, is unrealistic from many points of view. 
As to our solution (|165D , the second important aspect is that the parameters 
Xi{k) do not depend on frequency, since they depend essentially on elastic 
scattering processes. The dependence on frequency arises due to inelastic 
scattering processes which are contained in our self-energy operator. We 
proceed now with the derivation of the explicit expression for the self-energy. 
To calculate a high-order GF on the r.h.s. of (|158| ), we should use the second 
time variable (t') differentiation of it again. If one introduces the irreducible 
parts for the right-hand-side operators by analogy with the expression ( |160| ) , 
the equation of motion (|158|) can be rewritten exactly in the following form 



(174) FkA^) = fUu;) + FUu;)PkA^)FU^) 
Here the scattering operator P (36) is of the form 

(175) P,.(.;) = « ^tl\^tS 

Im 

In accordance with the definition (37), we write down the Dyson equation 

(176) F = F° + F°MF 

The self-energy operator M is defined by eq. (37). Let us note again that 
the self-energy corrections, according to (38), contribute to the full GF as 
additional terms. This is an essential advantage in comparison with the 
"Hubbard HI" solution and other two-pole solutions. It is clear from the 
form of Roth solution (|55| ) that it includes the elastic scattering corrections 
only and does not incorporate the damping effects and finite lifetimes. 
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For the full GF we find, using the formal solution of Dyson equation (pq), 
that it is equal to 

GkA^) = ( ^(a - nl,M++(^)) + -^{b - nZ^M^-{u;)) 

(177) [det((FO,M)-i-M,.H)]-i 

After some algebra, we can rewrite this expression in the following form 
which is essentially new and, in a certain sense, is the central result of the 
present theory: 

oj - {n+E_ + n-E+) - L 



fl78) G = 

{uj — Ej^ — n~Li){uj — E_ — n+L2) ~ n^n+L-^L^ 

where 

Li{k, uj) = Xi{k) 



L2{k,u) = \2{k) - ^M--ik,u;)- 



n 


—a 


n. 


—a 


n. 


—a 


n 


+ 


— cr 


n 




— cr 


n 


4- 
— cr 


n 


+ 


— cr 



L4{k,uj) = X4{h 

(179) L{k,u) = X{k) + nl^nI^(M++ + M— - M-+ - M+~) 

Thus, now we have to find explicit expressions for the elements of the self- 
energy matrix M. To this end, we should use the spectral theorem again to 
express the GF in terms of correlation functions 

(180) M^fiu) ~< > 

For the approximate calculation of the self-energy, we propose to use the 
following trial solution 

< Z)(-)t(t)I)(-) >«< aL(i)«^. X > 

+ < alnaitX-a >< n'^^^{t)aia > +/? < ^L,-, (*)««. >< aj,(tX-a > 

+af3 < bl^j_^{t)aia >< a]„{t)ba_^ > 

(181) +a(5 < blj^^{t)bii^„ >< a]^{t)ai^ > 
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It is quite natural to interpret the contributions into this expression in terms 
of scattering, resonance-broadening, and interference corrections of different 
types. For example, let us consider the simplest approximation. For this 



aim, we retain the first contribution in (181) 

[imiU= / ;(exp(/3u;0 + l 

/ + 00 
— exp(ia;'i)Af~^ ^ exp[-ik{Ri - Rj)]tiitm 



ijlm 



(182) J du;in{uJi) exp{iu;it)gmia{uJi) (^-■^ImK!^'^ {uJi - w')^ ■ 

Here K'^j^{uj) =<< nf_^\nj_^ »uj is the density-density GF. It is worthy 
to note that the mass operator (182) contains the term tutmj contrary to the 
expression ( |147| ) that contains the term {/^ . The pair of equations (|182| ) and 



( |176| ) is a self-consistent system of equations for the single-particle Green 
function. For a simple estimation, for the calculation of the self-energy 
( |182| ), it is possible to use any initial relevant approximation of the two-pole 
structure. As an example, we take the expression (|5^). We then obtain 



[IMI]«^«^|6(fc-g)|2Kf 



q 



n-rr 1 — n_ 



UJ — U — e{k — q)n^cr uJ — (-{k — q){i — n^^) 



In the same way, one can use, instead (|54D, another initial two-pole solution, 
e.g. the Roth solution (|55|), etc. 

On the basis of the self-energy operator ( |183| ) we can explicitly find the 
energy shift and damping due to inelastic scattering of quasi-particles. This 
is a great advantage of the present approach. 

In summary, in this Section, we obtained the most complete solution to the 
Hubbard model Hamiltonian in the strongly correlated case. It has correct 
functional structure, and, moreover, it represents correctly the effects of 
elastic and inelastic scattering in a systematic and convenient way. The mass 
operator contains all inelastic scattering terms including various scattering 
and resonance broadening terms in a systematic way. The obtained solution 
(|178| ) is valid for all band filling and for arbitrarily strong but finite strength 
of the Coulomb repulsion. Our solution contains no approximations except 
those contained in the final calculation of the mass operator. Therefore, we 
conclude that our solutions to the Hubbard model in the weakly correlated 



case ( |143| ) and in the strongly correlated case ( 178 ) describe most fully and 
self-consistently the correlation effects in the Hubbard model and give a 
unified interpolation description of the correlation problem. This result is 
to be contrasted with Hubbard , Roth and many other results in which this 
interpolation solution cannot be derived within the unified scheme. 
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It is clear from the present consideration that for the systematic construction 
of the advanced approximate solutions we need to calculate the collective 
correlation functions of the electron density and spin density and the density 
of doubles, but this problem must be considered separately. 

8.3 Correlations in Random Hubbard Model 

In this Section, we apply the IGF method to consider the electron-electron 
correlations in the presence of disorder to demonstrate the advantage of our 
approach. The treatment of the electron motion in substitutionally disor- 
dered Ar^Bi_x transition metal alloys is based upon a certain generalization 
of the Hubbard model, including random diagonal and off-diagonal elements 
caused by substitutional disorder in a binary alloy. The electron-electron in- 
teraction plays an important role for various aspects of behaviour in alloys, 
e.g. for the weak localization |9^. The approximation which is used widely 
for treating disordered alloys is the single-site Coherent Potential Approx- 
imation (CPA) [9^. The CPA has been refined and developed in many 
papers {e.g. ||97[| , |98f| ) and till now is the most popular approximation for 
the theoretical study of alloys. But the simultaneous effect of disorder and 
electron-electron inelastic scattering has been considered for some limited 
cases only and not within the self-consistent scheme. 

Let us consider the Hubbard model Hamiltonian (^9|) on a given configura- 
tion of an alloy {i') 

(184) H'^"^ = Hi"^ + H^"^ 
where 

ia ijcr 

(185) Hi""^ = ^J2 Ui^^ni.^ 

ia 

Contrary to the periodic model (|49|), the atomic level energy e^, the hop- 
ping integrals f^j^, as well as the intraatomic Coulomb repulsion are here 
random variables which take the values , f^^, and If^, respectively; the 
superscript z^(;u) refers to atomic species (i^, = A,B) located on site i(j). 
The nearest-neighbor hopping integrals were only included . 
To unify the IGF method and CPA into a completely self-consistent scheme 
let us consider the single-electron GF ( |153| ) Gija in the Wannier represen- 
tation for a given configuration (z^). The corresponding equation of motion 
is of the form (for brevity we omit the superscript (z^) where its presence is 
clear) 

(w - Ei) « aiala'j^ »uj= Sij + E**n << "'na\a]^ »uj 

n 

(186) +Ui « ni_^aia\a}-^ »^ 
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In the present paper, for brevity, we confine ourselves to the weak correla- 
tion and the diagonal disorder case. The generalization to the case of strong 
correlation or off-diagonal disorder is straightforward, but its lengthy con- 
sideration preclude us from discussing it this time. 

Using the definition (30), we define the IGF for a given (fixed) configuration 
of atoms in an alloy as follows 

(187) (^'■) « ni.,ai,\a]^ »= 
« ni-aaia\a^j^ » - < Ui-a >« aia\a^j^ » 

This time, contrary to ( [163D , because of lack of translational invariance we 
must take into account the site dependence of < rii-a >• Then we rewrite 
the equation of motion ( |186D in the following form 

^[(W -Ei-Ui < ni_„ >)6ij - tin] « ana\a]a »i^ = 
n 

(188) + t/,((*'^) 

In accordance with the general method of Section 3, we find then the Dyson 
equation for a given configuration {v) 

(189) G,,^(^) = Gl^iu) + J2 G^,^,{u)M^,Uuj)Gnja{u;) 

mn 

The GMF GF G^^^ and the self-energy operator M are defined as 

HimaG^j^{Lu) = 6ij 



m 



Pmna — ^mna ~l~ ^ ^ ^miaGij^Pjna 
ij 

Hima — Ui < ^^i— o- ^^^im tim 

(190) Pmnaiuj) = Um{^'''^ « nm-aama\nn-aai„ »^j'^)Un 

In order to calculate the self-energy operator M self-consistently, we have to 
express it approximately by the lower-order GFs. Employing the same pair 
approximation as ( |146D (now in the Wannier representation) and the same 
procedure of calculation, we arrive at the following expression for M for a 
given configuration (v) 

ImG^^l^^ {ui ) ImG^^l_„ (0^2 ) ImG'ii^l^ (cjs ) ; 
dujiduj2duj^ (1 — n{uji))n{uj2)n{(jj-i) 



(191) R 



u; + a;i-W2-W3 n{L02 + oo-i - u)i) 



As we mentioned previously, all the calculations just presented were made 
for a given configuration of atoms in an alloy. All the quantities in our 
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theory (G, G^, P, M) depend on the whole configuration of the alloy. To 
obtain a theory of a real macroscopic sample, we have to average over various 
configurations of atoms in the sample. The configurational averaging cannot 
be exactly made for a macroscopic sample. Hence we must resort to an 
additional approximation. It is obvious that the self-energy M is in turn a 
functional of G, namely M = M[G]. If the process of making configurational 
averaging is denoted by G, then we have 

G = G° + G^MG 

A few words are now appropriate for the description of general possibilities. 
The calculations of G^ can be performed with the help of any relevant avail- 
able scheme. In the present work, for the sake of simplicity, we choose the 
single-site CPA[96|, namely, we take 

(-.Qr,-. ^0 ,N _ ly-i sr^ exp(tA;(igm - Rn)) 

Here e(A;) = X]n=i ^n.o exp(iA;i?„), z is the number of nearest neighbors of 
the site 0, and the coherent potential is the solution of the CPA 

self-consistency equations. For the A^Bi-x, we have 

^^Lo) = xe\ + (1 - x)el - {^A - ^^F^i^, ^li^i - ^1; 
(193) F^u;,^n = G'imA^) 

Now, let us return to the calculation of the configurationally averaged total 
GF G. To perform the remaining averaging in the Dyson equation, we use 
the approximation 

GHIG w G^MG 

The calculation of M requires further averaging of the product of matrices. 
We again use the prescription of the factorizability there, namely 



M PS (UmUn) (ImG) (ImG) (ImG) 



However, the quantities UmUn entering into M are averaged here according 
to 



UmUn = U2 + {Ui - U2)5mn 

Ui = x^Ul + 2x(l - x)UaUb + (1 - x fU^ 

(194) U2 = xUl + il-x)U^ 

The averaged value for the self-energy is 

(195) M-ranaiuj) 
U2 f 

= \ R{^^l,^^2,(^3)ImGnm-a{i^l)IrnGrnn-a{i^2)I'mGrnna{i^3) + 

^-^-l^^^^^mn J ^(Wl, W2, W3)/mG,im-CT(wi)/m(5mn-o-(w2)-/'"lG'mna(w3) 
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The averaged quantities are periodic, so we can introduce the Fourier trans- 
form of them, i.e. 

k 

and similar formulae for G and . Performing the configurational averaging 
of the Dyson equation and Fourier transforming of the resulting expressions 
according to the above rules, we obtain 

(196) Gkaiuj) = [10- e{k) - - Mk.{u)]-^ 
where 

-MfecrCo;) = V / R{LOl,i02,i^z)N~'^ImGp^g-.„{uJl)ImGq^„{L02) 

(197) \U2lmGk+pAu:s) + ^ ImGk+p-g{oo3)] 

9 

The simplest way to obtain an explicit solution for the self-energyM is to 
start with a suitable initial trial solution as it was done for the periodic 
case. For a disordered system, it is reasonable to use, as the first iteration 
approximation the so-called Virtual Crystal Approximation (VGA): 

— /mGL^^(u; + ie)p.5(u;-i?,-) 
vr 

where for the binary alloy A^Bi^x this approximation reads 
y = xV"^ + {I - x)V^ ; = + e{k)- 



Note, that the use of VGA here is by no means a solution of the correlation 
problem in VGA. It is only the use of the VGA for the parametrization of 
the problem, to start with VGA input parameters. After the integration of 

(197) the final result for the self-energy is 

(198) Mfca(w) = 
U2 ^ n{E^^^)[l - n{E--) - n{El^^)] + n{El^^)n{E-") 



+ 



h ^ + ^pT, - Er - El^p_g 



It is to be emphasized that the equations (|195| ) - ( |198| ) give the general mi- 
croscopic self-consistent description of inelastic electron-electron scattering 
in an alloy in the spirit of the GPA. We took into account the randomness not 
only through the parameters of the Hamiltonian but also in a self-consistent 
way through the configurational dependence of the self-energy operator. 
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8.4 Electron-Lattice Interaction and MTBA 



To understand quantitatively the electrical, thermal, and superconducting 
properties of metals and their alloys, one needs a proper description of an 



electron- lattice interaction too |9^, 1 102], [103]- A systematic, self-consistent 



simultaneous treatment of the electron-electron and electron-phonon interac- 
tion plays an important role in recent studies of strongly correlated systems. 
It was argued from different points of view that to understand quantitatively 
the phenomenon of high-temperature superconductivity one needs a proper 
inclusion of electron-phonon interaction, too. A lot of theoretical searches 
for the relevant mechanism of high temperature superconductivity deal with 
strong electron-phonon interaction models. The natural approach to the 
description of superconductivity in that type of compounds is the modi- 



fied tight-binding approximation (MTBA) |9g], []lO||. The papers 0,|9|] 



[|100(| ,|101] contain a self-consistent microscopic theory of the normal and su- 
perconducting properties of transition metals and strongly disordered binary 
alloys in the framework of the Hubbard Model (|^) and random Hubbard 



model (184). Here we derive a system of equations for the superconductiv- 
ity for tight-binding electrons of a transition metal interacting with phonons 
within the IGF approach. We write the total Hamiltonian of the electron-ion 



system as the sum |71] 



(199) H = He + Hi + He-i 

where Hg is the electron part of the Hamiltonian represented by the Hubbard 
operator (^9|). The Hamiltonian of an ion subsystem and the operator of 
electron-ion interaction have the form 

(200) h, = ^Y.§^ + IT. Ki<ui 

^ mnaf) 

(201) H,., = E E ^^"(^^>l«.-< 
where 

(202) J2 y^jiRnX = - n,) 

Here Pn is the momentum operator, M is the mass of an ion, and Un is the 
displacement of the ion from the equilibrium position at the lattice site Rn- 
In a more convenient notation the electron-phonon interaction Hamiltonian 
in the modified tight-binding approximation reads 

(203) He-i = E E ^"(^' ^ + Q)QgA+qaaka 
where 

(204) (fc, + = jj^^^ E (^) [sin Sj- sin a„ (k - q}] 
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here qo is the Slater coefficient |^9| having the origin in the exponential 
decrease of the wave functions of d-electrons, N is the number of unit cells 
in the crystal, and M is the ion mass. The quantities e,^{q) are polarization 
vectors of the phonon modes. 
For the ion subsystem, we have 

(205) Hi = l ^(P|,P,. + u;\qi.)QlQ,,) 

where Pqy and Qqy are normal coordinates, and u;(gi') are acoustical phonon 
frequencies. It is important to note that in spite of the fact that in Hubbard 
model (^9|) , the d- and s(p)-bands are replaced by one effective band , the s- 
electrons give rise to screening effects and were taken into effects by choosing 
the proper values of \J and the acoustical phonon frequencies. 

8.5 Equations of Superconductivity 

To derive the superconductivity equations, we use the IGF method of Section 
3 in which the decoupling procedure is carried out only for approximate 
calculation of the mass operator of the matrix electron GF. According to 
the arguments of Section 4.3 , eqn. (|6^) , the relevant matrix GF is of the 
form 

As was discussed in Section 4.4, with this definition, one introduces the so- 
called anomalous (off-diagonal) GFs which fix the relevant BCS-Bogoliubov 
vacuum and select proper symmetry broken solutions. Differentiation of 
Gij{t — t') with respect to the first time gives for the Fourier components of 
the equations of motion 

(207) ^{^^ij - tij) « aja\4'a »= ^ii' + 

j 

U « aiuni-a\a\,^ » +^Kin << ajaUn\a\,^ » 

nj 

(208) ^(cjJi,- + %) « a]_Jal^ »= 

j 

-U « a\_„ni„\a\,^ » +^ V,m « «]-a^n|oL » 

nj 

Following the general strategy of the IGF method, we separate the renor- 
malization of the electron energy in the Hartree-Fock-Bogoliubov generalized 
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mean field approximation (including anomalous averages) from the renor- 
malization of higher-order due to inelastic scattering. For this, we introduce 
irreducible parts of the GF in accordance with the definition ( as an example, 
we take two of the four Green functions) 

(209) ((^'')<< ai^a\_^ai^„\a\ ^ >>^) =<< ai^a\_^ai^^\a\ ^ »^ - 

- < Ui^a- > ^11+ < Qiatti-a >« al_Jal,^ »^ 
(^^'^« 4aaia4^Jal^ »^) =« alai^aj^Jal^ »^ - 

From this definition it follows that this way of introducing the IGF broadens 
the initial algebra of the operators and the initial set of the GFs. This 
means that "actual" algebra of the operators must include the anomalous 
terms from the beg inning, naniGly: (tt^cr? ^lo-' ^icrt ^Jfj^l— tr' ^i—o^ia 

). The 

corresponding initial GF is the form ( p06D . The choice of the irreducible 
parts of the GF in ( ^091) is specified by the "orthogonality" constraint (31), 
which makes it possible to introduce unambiguously the irreducible parts 
and make the inhomogeneous terms in the equations for them vanish. Using 
( |209| ) , we rewrite eqs.( |207| ) and ( |208| ) in the form 

(210) ^({uj-U < >)6,j - tij^ « ajvl4<T >>= ^n' 

3 

(211) ^((w + [/ < Uj^ >)6ij + tji^ « 4_Jal^ »= 

3 

+u < aL«,U >« «*-i4<x » "E « (pji— 4— )^"'^i4a » 

3 

where 

(212) pij„ = Unj_„5ij + VijnUn{l - %) 



In the representation of the Nambu operators |71) 



(213) ^i,-. = \y] = 



the equation of motion for GF ( ^lOf ) can be represented as 

(214) ^(wroJij - tijT3 - S^^) « %bj\ipl »= 

3 

Sii'To + J2 « (ftj7-3V'i)^^''^|V'l' » 
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Here the Hartree-Fock-BogoHubov elastic Coulomb term ( |6^ is of the form 
(215) Ef, = -Ut3 < > ^3 + |(ro + T3) 



To calculate the irreducible matrix GF in ( |214| ), we write down for it the 
equation of motion with respect to the second time t' and then separate the 
irreducible part with respect to the operators on the right-hand-side of the 
corresponding GF. This gives the Dyson equation in the matrix form 

(216) Gu'iu:) = G%{u:) + G%{u:)M,j,{u:)Gj,,{u:) 

jj' 

The generalized mean field GF G^ and the mass operator are defined by 

(217) ^(wTo<5,j - UjTs - S^^)G°i, = 6,i'To 

j 

(218) Mkk' = E(« iPk,r3i^jf'^\{ij],T3p,>k'y^'''^ »)(?) 

jj' 

The explicit expression for the mass operator ( |218D is of the form 

(219) Mii,{u;) = 

((-)« a,tPi,T|p,Y^aJ,^ »(-))b) ((-)« a,^p,,^\pj>,>^ay^ »(-))(p) ' 



E 



n 



The mass operator ( |2ig| ) describes inelastic scattering of electrons ( the 
elastic part is contained in T,'^^ ) on fluctuations of the density of a total 
electron- ion charge in the lattice. To find an approximating expression for 
the mass operator ( pl9D , we adopt the following trial approximation 



(220) < Pj'i'cr{t)a'-,^{t)ajcrPija >^"'^W< Pj>i'a{t)pija >< <4'a(^)"'j<^ > 

This approximation was made in the spirit of the approximation of "two in- 
teracting modes" and means ignoring the renormalization of the vertex, i.e., 
the correlation in the propagation of an electron (hole) and the propagation 
of charge density fluctuations. 

Writing down further spectral representation for the correlation functions in 
( |220| ) , we represent the mass operator by the sum 



(221) M,i, {u) = Ml, ioj) + Ml {uj) 

The first contribution has a form characteristic of an interacting electron- 
phonon system 

^iA^) = 1^1^ VijnVj'i'n'i: / (cot — + tan — ) 

nn' jf 

(222) {--Im « Un\Un' »^,){--T3lm « ^^^#1 »uji T3) 



TT vr 
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The contribution Mf^, has a more comphcated structure 

(223)M^, = ^P ^^^^"^ (cot^+tan^)f"^^^ "^^^ 
2 J-oo u; - - W2 2 2 ' V m2i ^22 

where 

mil = {-—Im « nil |nj/| >>tj2)( -^"^ << aitl'^I'T >>'^i) 

vr vr ' 

W.12 = {-Im « nji|nj/t >>a;2)(--/m << ajt|aj,| >>(^i) 
vr vr ■'' 

m22 = {--Im « nj||nj/t »^,^){--Im « aii\ai: 
vr vr 



The equations (|216| ) and (|221|) constitute a self-consistent system of equa- 
tions for the single-particle GF of the Hubbard model on a vibrating lattice. 
Note that these equations of superconductivity can be in an obvious way 
transformed to the standard form of the Eliashberg equations [^]. The 
numerical calculations of the electron-phonon spectral function a'^{uj)F{Lo) 
for a few transition metals were done in ref. |102[|. It is worthy to empha- 



size that in paper [101] a very detailed microscopic theory of the strong 
coupling superconductivity in highly disordered transition metal alloys was 
developed on the basis of the IGF method within the MTBA reformulated 
approach [|10[)(| . The Eliashberg-type strong coupling equations for highly 
disordered alloys were derived. It was shown that the electron-phonon spec- 
tral function in alloys is modified strongly. Thus, the self-consistent sys- 
tem of superconductivity equations obtained in the Wannier representation 
makes it possible to investigate real transition metals, their alloys, and com- 
pounds from a unified point of view. 



9 Quasi-Particle Dynamics of Anderson Models 
9.1 Quasi-Particle Dynamics of SIAM 

In this Section, we consider the many-body quasi-particle dynamics of the 
Anderson impurity model at finite temperatures in the framework of the 
equation-of-motion method. In spite of many theoretical efforts, there is no 
complete solution of the dynamic problem for the "simple" Anderson/Hubbard 
model. One of the main reasons for this is that it has been recognized rel- 
atively recently only that the simplicity of the Anderson model manifests 
itself not in the many-body dynamics ( the right definition of quasi-particles 
via the poles of GF; see Section 6.1 ) but rather at quite a different level 
- in the dynamics of two-particle scattering, resulting in the elegant Bethe- 
ansatz solution (for the relativistic spectrum linear in k ), which gives the 
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static characteristics (static susceptibility, specific heat, etc). In this sense, 
as to the true many-body dynamics, the complete analytic solution of this 
problem is still quite an open subject. This Section is primarily devoted 
to the analysis of the relevant many-body dynamic solution of the SIAM 
and its correct functional structure. We wish to determine which solution 
actually arises both from the self-consistent many-body approach and in- 
trinsic nature of the model itself. We believe highly that before numerical 
calculations of the spectral intensity of the Green function at low energy 
and low temperature it is quite important to have a consistent and closed 
analytic representation for the one-particle GF of the SIAM and Hubbard 
model. The paper [|104|| clearly shows the importance of the calculation of 
the GF and spectral densities for SIAM in a self-consistent way. An alterna- 
tive approach to dynamics of the Anderson model was formulated within a 
modified version of the Kadanoff-Baym method [105], |10(:]. Unfortunately, 
the Neal approach also have certain drawbacks. 

A proper many-body description of dynamic correlations is very actual also 
for the investigation of the dynamics of the many-impurity Anderson model, 
where standard advanced many-body methods do not work properly in usual 
formulation. Recently, a lot of efforts were devoted to a better understanding 
of the static and dynamic properties of the Anderson model in the context of 
many- impurity case ||2^. This field is quite important for the description of 
magnetic properties of anomalous rare-earth compounds [^],||6^. The prob- 
lem of an adequate and consistent description of dynamics of single-impurity 
and many-impurity Anderson models ( SIAM and MIAM ) and other mod- 
els of correlated lattice electrons was not yet solved analytically completely 
. During the last decades, a lot of theoretical papers were published, at- 
tacking the Anderson model by many refined many-body analytic methods. 
Nevertheless, a fully consistent dynamic analytic solution in the closed form 
for a single-particle propagator of SIAM is still lacking. In this Section, the 
problem of consistent analytic description of the many-body dynamics of 
SIAM is discussed in the framework of the equation-of-motion approach for 
double-time thermodynamic GFs. In addition to the IGF approach, we find 
a new exact identity relating the one-particle and many-particle GFs. Using 
this identity, we present a consistent and general scheme for construction 
of generalized solutions of SIAM. A new approach for the complex expan- 
sion for the single-particle propagator in terms of Coulomb repulsion U and 
hybridization V is proposed. Using the exact identity, an essentially new 
many-body dynamic solution of SIAM is derived. This approach offers a new 
way for the systematic construction of approximate interpolation dynamic 
solutions of strongly correlated electron systems. 
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9.2 IGF Approach to SIAM 

After discussing some of the basic facts about the correct functional struc- 
ture of the relevant dynamic solution of correlated electron models we are 
looking for, described in previous Sections, we give a similar consideration 
for SIAM. It was shown in |25], using the minimal algebra of relevant opera- 
tors, that the construction of the GMFs for SIAM is quite nontrivial for the 
strongly correlated case, and it is rather difficult to get it from an intuitive 
physical point of view. Let us consider first the following matrix GF 



(224) G{uj) 



« /0a|4<x » « /oal/oV » J 

Performing the first-time differentiation and defining the irreducible GF 

(225) ((-)« foafljo-a\fl =« f0aflj0-a\fl ». " 

- < no^a >« foalfoa »uj 

we obtain the following equation of motion in the matrix form 

(226) J2 Fp{^)Gp{oo) = i + UD^"'\lo) 

p 

where all definitions are rather evident. Proceeding further with the IGF 
technique, the equation of motion (226) is exactly rewritten in the form of 
the Dyson equation 

(227) G{uj) = G°(w) + G^{uj)M{uj)G{uj) 
The generalized mean field GF G^ is defined by 

(228) Y.WGl{u;)=i 

p 

The explicit solutions for diagonal elements of G^ are 

(229) « Mfl »»= (u, - E|,„ - C/n_„ - S{ui))y' 

(230) «c.„|cL»»=(.-.»-;J!!^)- 
where 

(231) ^(-^'"^i^ 

k 

The mass operator, which describes inelastic scattering processes, has the 
following matrix form 





(232) M{cv) 



Mo, 
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where 

(233) Mo. = fo.no-MLno-a 

From the formal solution of the Dyson equation (38) one obtains 

(234) « fo^lfl »^= {lj - Eoa - Un.„ - Mq, - 5(l^))"' 



(235) << Cfc.|4^ >>(^= (w - efc 



to - E'oa - U - Mq^ 



To calculate the self-energy in a self-consistent way, we have to approximate 
it by lower-order GFs. Let us start by analogy with the Hubbard model 
with a pair-type approximation ( |146D 



(236) Mo^iuj) 

dujiduj2duj'i 



C/2 



OJ + UJl — UJ2 — OJ^ 

where we used the notation 

9o<7(w) = --Im « foalfL »^ 

The equations ( |227|) and ( ^36|) constitute a closed self-consistent system of 
equations for the single-electron GF for SIAM model, but only for weakly 
correlated limit. In principle, we can use, on the r.h.s. of ( ^361) , any workable 
first iteration-step form of the GF and find a solution by repeated iteration. 
If we take for the first iteration step the expression 



(237) 50a (w) ~ 5(w - Eq, - Un. 

we get, for the self-energy, the explicit expression 

Mo.(.) = ^2 n(ii;o. + ^n„.)(l-yo. + ^n_.)) ^ u'N.M-N-.)Gliu.) 
(238) 

where A^-o- = ti(E()„ + Un^a)- This is the well-known "atomic" limit of the 
self-energy. 

Let us try again another type of the approximation for M. The approxi- 
mation which we will use reflects the interference between the one-particle 
branch and the collective one 

< h.{t)flAt)h~a{t)fljo-.fl >(^^)« 

< f Lit) ha >< no-a(t)?^0-a > + 

(239) < fl_At)k-a X fo-At)foAt)fLfLa > 
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If we retain only the first term in ( p39D and make use of the same iteration 
as in (12371), we obtain 



(240) Mo.(.) - ^2 (l-n(^o. + ^n..)) ^ ^ 



If we retain the second term in ( 239 ), we obtain 



. X 9 /■+°° , 1 + N(uJi) - n(uJ2) 



oo 



u; — wi — c<j2 



(-i/m << S^\S^ 
vr 

(241) (--Im « fo.\fl », 

^ vr 

where the foUowing notation was been used: 

Sq = /d|/oi; Sq = /o|/oT 
It is possible to rewrite (241) in a more convenient way now 



Afo.(w) = y do;' (cot + tan ^) 

(242) (_limxT±(^ _ uj')go^{Lo')) 



The equations ( |227| ) and ( |242| ) constitute a self-consistent system of equa- 
tions for the single-particle GF of SIAM. Note that spin-up and spin-down 
electrons are correlated when they occupy the impurity level. So, this really 
improves the H-F theory in which just these correlations were missed. The 
role of electron-electron correlation becomes much more crucial for the case 
of strong correlation. 

9.3 SIAM. Strong Correlation 

The simplest relevant algebra of the operators used for the description of 
the strong correlation has a similar form as for that of the Hubbard model 
( |154| ). Let us represent the matrix GF (|224| ) in the following form 



(243) Giu) = Y: ( 

al3 \ 



« Cka\c\^ » « Cka\dli3^ » 
« doaa\4-^ » « doaa\dL^ » 



Then we proceed by analogy with the calculations for the Hubbard model. 
The equation of motion for the auxiliary matrix GF 

(244) F^ico) = 

f « Cka\c\„ » « Cfc<xl4+a >> « Cfc<7l4-a >> 

<< do+(7l4<7 >> << dQ+Mlj^^ » « (io+<x|4-(7 >> 
V<< do-a|4<T >> <<do-al4+a>> «do-a\dl_^». 

75 



is of the following form 

(245) EF^{uj) -i = D 
where the following matrix notation was used 

(246) E = 





-Vk 


-Vk 





(a; - Eo, - U+) 











{uj - Eo^ - [/_) 

no 








VO no_^ 



U, a = + 
0, a = — 



Here Z) is a higher-order GF, with the following structure 

/ 

(247) D{lo) = D21 D22 D23 

\D3l Di2 D33. 

In accordance with the general method of Section 3, we by define the matrix 
IGF: 

(248) D(^^Hu;) = D-Y.(^f_iyG:^ G^) 
Here the notation was used: 

(249) A++ = ^ (/o-g^P"'^ + (i-afo-<T)inoa - no_^) > 

< no-<T > 

f250l A— = ~ ^ ^■^o-'^^P-'" ^ 4-'t/o-'^)(1 + np^ - np,^) > 

< 1 - rao-<7 > 

A~+ = A++, A+- = -A— 
The generalized mean-field GF is defined by 

(251) ^F°(a;)-i = 0; = J] 

From the last definition we find that 

f I ft ^^o_ < "0-^ > n I T.pVp^"^ ^ 

« ^_E,^-U+- Y., V,A++ ^^^u-Eo^-uJ 

(252) l-<no_^> . Ep'^p-^^' ^ 
^ ' ^u-Ep^-U.-Y^pVpA—^ ^ uj-Eq^-uJ 

(253) « Cka\c\, »l= {oJ-ek- \Vk\'' F^^^uj))'^ 
where 



(254) F"* = _1^:}!ZIL^ + 

W - E'oa -U+ UJ - Eoa - U- 

For 1^ = 0, we obtain, from solution ( p52[) , the atomic solution F'^^. The 
conduction electron GF (253) also gives a correct expression for Vk = 0. 
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9.4 IGF Method and Interpolation Solution 

To show explicitly the flexibility of the IGF method, we consider a more 
extended new algebra of operators from which the relevant matrix GF should 
be constructed to make the connection with the interpolation solution of 
the Anderson model. For this aim, let us consider the following equation of 
motion in the matrix form 



(255) 



where G is the initial 4x4 matrix GF and D is the higher-order GF: 



(256) 



Gcr 



( G\\ G\2 G\z 

G21 G22 G23 

Gsi G32 G33 

\ G41 G42 G43 



Gl4 \ 

G'24 

^34 
G44 / 



Here the following notation was used 



G 



11 



= << Cka\cl^ >>; G12 =« CkalfL » 



(257) 



Gi3 =<< Cfc^l/oV^-o-fT >>; Gu =« Ckalc'k^no^a » 

G2I =« /0<xl4a »; G22 =« foalfL » 
G23 =<< focrlfoa^O-cT >>; G24 =« /o<7|4ct^O-(7 >> 
G3I =<< foanO-a\cl^ >>; G32=« foanO-a\foa » 
G33=« foano-a\foa'^0-a >>; G34 =<< /oa^^O-CT | C^a'^O-a >> 
G41 =<< Cfccrno-cr|4<7 '^42 =<< CkanO-a\fL » 

43 =<< Ckano-alfo^no^a- »] G44, =« CkanQ-a\c\^nQ^a^ »] 



G. 



We avoid to write down explicitly the relevant 16 GFs, of which the matrix 
GF D consist, for the brevity. For our aims, it is enough to proceed forth 
in the following way. 

The equation (255) results from the first-time differentiation of the GF G 
and is a starting point for the IGF approach. Let us introduce the irreducible 
part for the higher-order GF D in the following way 



(258) 



= - E L^'^Gafs; («, /3) = (1, 2, 3, 4) 



and define the GMF GF according to 



(259) 



J2F(j^,k)G^'/{u;) = I, 
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Then, we are able to write down exphcitly the Dyson equation (37) and the 
exact expression for the self-energy M in the matrix form: 



(260) 


























M33 


M34 


Vo 





M43 


M44/ 



r-1 



Here the matrix I is given by 







1 

1 < no_a > 

< no-CT > < no-a > 
V < no-<7 > 

and the the matrix elements of M are of the form: 

M33 = (« A^rHpM'^'Hq) »)(^),M34 = (« 4^\p)\Bt\k,q) »] 








M43 = (« Af\k,p)\Br\q) »)(^),M44 = (« Af\k,p)\Bi'^\k,q) » 



ip) 
\ip) 



(261) 



where 



(262) 



Blip) = ifL4y-ak-a - fljl_^C^-a 



t A 



B2ik,p) = (4aCj-<T/0-a - C^/o-aCp-a); 



t #t 



Since the self-energy M describes the processes of inelastic scattering of 
electrons (c-c , f-f, and c-f types), its approximate representation would be 
defined by the nature of physical assumptions about this scattering. 
To get an idea about the functional structure of our GMF solution ( ^59|) , 
let us write down the matrix element G^^^: 



^33^ =« kano-a\foa'^0-a » = 



UJ 



U-S^'^P{uj)-Y{io] 



+ 



(263) 
(264) 



[UJ 



MF 



-U -S^'^P{io)-Y{Lo)){io-Eoa-S{u;)) 

UZ{uj) 



Z{co) = S{io)Y^ 



UJ 



cMF 



+ 



(265) 



UJ 



MF 



32 
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Here the coefficients L"^^ , L^^ , L^^ , and are certain complicated averages 
(see definition (|25§| )) from which the functional of the GMF is build. To 
clarify the functional structure of the obtained solution, let us consider our 



first equation of motion ( 255 ) , before introducing the irreducible GFs (^55 



Let us put simply, in this equation, the higher-order GF D = 01 To dis- 



tinguish this simplest equation from the GMF one ( 259 ), we write it in the 
following form 

(266) ^F(p,fc)GO(p,ai) = / 

p 

The corresponding matrix elements which we are interested in here read 

(267) G°2 =« fo^lfl »= 

1- < no-(T > , < no-a > 



uj — Eq„ — S{uj) uj — Eocr — S{uj) — U 



(268) =« /o.no-.|/t.no-. »= - _ _ ^^^^ _ ^ 

(269) =« foano-a\fL »= Gl, 

The conclusion is rather evident. The simplest interpolation solution follows 
from our matrix GF ( p56| ) in the lowest order in V, even before introduction 
of GMF corrections, not speaking about the self-energy corrections. The 
two GFs G32 and G33 are equal only in the lowest order in V. It is quite 
clear that our full solution (38) that includes the self-energy corrections is 
much more richer. 



It is worthwhile to stress that our 4x4 matrix GMF GF (256) gives only ap- 
proximate description of suitable mean fields. If we consider more extended 
algebra, we get the more correct structure of the relevant GMF. 



9.5 Dynamic Properties of SIAM 

To demonstrate clearly the advantages of the IGF method for SIAM, it 
is worthwhile to emphasize a few important points about the approach 
based on the equations-of-motion for the GFs. To give a more instruc- 
tive discussion, let us consider the single-particle GF of localized electrons 
Ga =<< /oo-l/oo- >>• The simplest approximate "interpolation" solution of 
SIAM is of the form: 

Ga{0J) = ^r—- + 



UJ - Eo„ - S{oj) {uj - Eq„ - S{uj) - U){uj - Eo„ - S{uj)) 
(270) l-<m.^> ^ < no-. > 



UJ — Eocr — S{uj) UJ — Eocr — S{uj) — U 
The values of are determined through the self-consistency equation 

(271) =< no. >= -- [ dEn{E)ImG„{E,n„) 

vr J 
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The " atomic- like" interpolation solution ( pTOD reproduces correctly the two 
limits: 

G (w) = 1~ < '^o-^ > j_ < ^o-<T > V = {) 

^ — Eqo ^ — Eq(j — U 

(272) G.H = p ^ /or = 



The important point about formulas ( |272| ) is that any approximate solution 
of SIAM should be consistent with it. Let us remind how to get solution 
( |272| ). It follows from the system of equations for small-V limit: 

(W - Eoa - Siuj)) « foalfL »^= ^ + U « /oano-al/oV 

{UJ - Eoa -U) « foano-alfoa >>a;~ 

(273) < no-a > + ^ « CkanOalfL »UJ, 

k 

{uJ - ek) « CkanO-a\foa »lo = 
Vk « foanO-a\fL >>'^ 

The equations ( |273D are approximate; they include two more terms, treated 
in the limit of small V in paper |107|]. 



We now proceed further. In paper [107] the GF G was calculated in the 
limit of infinitely strong Coulomb correlation U and for small hybridization 
V. The functional structure of the Lacroix solution generalizes the solution 
(|272| ). The starting point is the system of equations: 



(274) {uj - Eo^ - S{lo)) « foML»=^ + U « hano-^\fl » 
(a; - Eo^ - U) « foano-a\fo^ »=< no-a > Ckano-alfoa » - 

k 

(275) « Ck-jl_jQa\fL » + « 4-a/o<x/o-a|/oV » 

Using a relatively simple decoupling procedure for a higher-order equation 
of motion, a qualitatively correct low-temperature spectral intensity was 
calculated. The final expression for G for finite U is of the form 

« foAfl »= ^_E^^-S{u) + USi{u:) + 
^27g) C/<no-. >+C/Fi(a;) 



K{iv){LO - Eo^ - S{u;) + USiiio)) 

where Fi, Si, and K are certain complicated expressions. We write down 
explicitly the infinite U approximate GF |107|: 

(2T7) « MfL 1- < "0-^ > 



u) - Eo„ - S{uj) - Zl{u)) 
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The following notation was used 

(278) F^ = vy < fo-aCk-a > 

k --^^ 

t t 

(279) = y ^ C^-aCfc-a > _ g^^^yy < fo-a^k-a > 

t^k ^-'^ k ^-'^ 

The functional structure of the single-particle GF ( |276| ) is quite transparent. 
The expression in the numerator of (p76D plays the role of "dynamic mean 



field", proportional to < fl^^Ck-a >■ In the denominator, instead of bare 



shift S{uj) ([23TD we have an "effective shift" = S{oj) + Zl{uj). The choice 



of the specific procedure of decoupling for the higher-order equation of mo- 
tion specifies the selected "generalized mean fields" (GMFs) and "effective 
shifts". 

9.6 Interpolation Solutions of Correlated Models 

It is to the point to discuss briefiy the general concepts of construction of an 
interpolation dynamic solution of the strongly correlated electron models. 
The very problem of the consistent interpolation solutions of the many-body 
electron models was formulated explicitly by Hubbard in the context of the 
Hubbard model. Hubbard clearly pointed out one particular feature of con- 
sistent theory, insisting that it should give exact results in the two opposite 
limits of very wide and very narrow bands. The functional structure of a 
required interpolation solution can be clarified if one considers the atomic 
(very narrow band) solution of the Hubbard model (p9|): 

(280) G^'ioj) = + ''-'^ ^ 



w — to uj — to — U uj — to — S"*(a;) 
where 

(281) ^"M= ^ g'-f.,) !; ; to = 

Ul-to 

Let us consider the expansion in terms of U: 

(282) S'**(aj) w n_^C/ + n_^(l - n^„)U^^— + 0{U) 

UJ - to 

The "Hubbard I" solution (|5^can be written as 

1 1 



(283) Gfc 



Lu - e{k) - S'^*(w) (C^*)-! + to- e{k) 



The partial "Hubbard HI" solution, called the "alloy analogy" approxima- 
tion is of the form: 

(284) S(a;) ""'"^ 



1 - ([/ - S(w))G(w) 
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Equation (284) follows from (pSlj) when one takes into account the following 



relationship: 

(285) oc G(uj) - S(w)G(aj) 

a; — to 1 — ^-o- 

The Coherent Potential Approximation (CPA) provides the basis for physi- 
cal interpretation of equation ( p84D which corresponds to elimination of the 
dynamics of —a electrons. By analogy with ( |282D , it is possible to expand: 

(286) frr'^Z]^ , « ^--^ + n_^U{U - ^)G\u; - S) + 0{U) 
1 — [U — L,{uj))G{uj) 



The solution ( PTTf ) does not reproduce correctly the U-perturbation expan- 
sion for the self-energy: 



J dEi j dE2 j dEs 



M„{uj) -^U < no-a > + 
n{Ei)n{E2){l - n{E^)) + (1 - n(^i)(l - n{E2))n{E^) 
uj — El — E2 + E^ 
(287) ImG„{Ei)ImG-„{E2)ImG-„{E^) 

It can be shown that it is possible , in principle, to find a certain way to incor- 
porate this U"^ perturbation theory expansion into the functional structure 
of an interpolation dynamic solution of SIAM in a self-consistent way within 



the higher-order GFs [108|. The IGF approach with the use of minimal al- 
gebra of relevant operators allows one to find an interpolation solution for 
weak and strong Coulomb interaction U and to calculate explicitly the quasi- 
particle spectra and their damping for both the limits. The U-perturbation 
expansion ( |147D is included into the IGF scheme in a self-consistent way. 
The correct second-order contribution to the local approximation for the 
Hubbard model is of the form 

(288) oc r^o-An, » 

n_^(l - n_^j 

The same arguments are also valid for SIAM. 

9.7 Complex Expansion for a Propagator 

We now proceed with analytic many-body consideration. One can attempt 
to consider a suitable solution for the SIAM starting from the following exact 
relation derived in paper |28]: 

(289) «foa\fl»=g' + g'P9° 

(290) 5° = (c^ - Eo^ - S{u;))-^ 

(291) P = U < no-a > « focrnO-alfL^O-a » 



The advantage of the equation ( 289D is that it is a pure identity and does 



not include any approximation. If we insert our GMF solution (277) into 
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( |289| ) , we get an essentially new dynamic solution of SIAM constructed on 
the basis of the complex (combined) expansion of the propagator in both 
U and V parameters and reproducing exact solutions of SIAM for V = 
and U = 0. It generalizes (even on the mean-field level) the solutions of 
papers gO^, p5| . 

Having emphasized the importance of the role of equation (|289|) , let us see 



now what is the best possible fit for higher-order GF in ( 2911 ). We consider 
the equation of motion for it: 



(292) 



(CJ - -U) « foano~a\fL'>T'0-a » = < TlQ-a > + 

^ Vfc(<< Cfea-^O-al/oa^O-a >> + 



« cl_^foafo^a\foa^O-a » " << Ck-afo-afoa\fL'^0-a >>) 

We can think of it as defining new kinds of elastic and inelastic scattering 
processes that contribute to the formation of generalized mean fields and 
self-energy (damping) corrections. The construction of suitable mean fields 
can be quite nontrivial, and to describe these contributions self-consistently, 
let us consider the equations of motion for higher-order GFs in the r.h.s. of 
(|293) 



{to - ek) « Cfccr'T-o-al/oa^o-fj >>= 

(293) V « foam-a I fL^O-a » + 



(294) 



(a; - efc - Eq^ + £^0-a) << Ck-afo-afoa\fL'>T'0-a » 

= - < fo-a'^k-anoa > " 
V « foano~a\foanO~a » + 
^y(<< Ck-afl^a'^pAfL^O-a » - « Ck^acl^af0a\fL^T'0-a >>) 



(295) {uj + ek- Eq„ - Eo-<7 -U)« cl_^hak-a\fLnQ^^ » 

= - < > + 

V « fo^no^a\foa^O-a » + 
Ck-a^P'^fo-alfLnO-a » + « cl_^foaCp-a\foanO-a >>) 

P 

Now let us see how to proceed further to get a suitable functional structure 
of the relevant solution. The intrinsic nature of the system of the equations 
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of motion ( |293D - ( |295| ) suggests to consider the following approximation: 

(296) (CJ -ek) « CkanO-a\foa'^0-a »^V « foanO-a\foa'^0-a » 
- efc - Eoa + Eo^^) « Ck-afo-afoa\foa^O~a >>» " < fo-aCk-anQa > 

(297) -V{« /o<xno-<x|/o<,rao„a » - « Cfc_^4_^/o<x|/oV^o-a ») 

{uj + ek- Eoa - Eo-a -U) « 4_^/oa/o-a|/oa^O-<7 »^ - < cl_^foafoafo-a > + 

(298) V{« /oano-<x|/o<,no_^ » + « cl_^foaCk-a\fL^O-a ») 

It is transparent that the construction of approximations ( |296| ) - ( |298| ) is 
related with the small-V expansion and is not unique, but very natural. As 
a result, we find the explicit expression for GF in (291) 



(299) « foano^alfL^O-a » 

Here the following notation was used 
(300) 

k 

(301) 



UJ 



1 

efc — -E'oo- + E{ 



+ 



O-cr 



u - Eq^ - U - Si{uj) 

Si{u)) = S{uj) 

) 

OJ + Ck — Eqct — Eo^cr — U 

Fl = Y.{VF2 + V^F^) 

k 
ft 



(302) 
(303) 



< cl_^foafoafo-a > ^ < f^_„Ck-anOa > 



UJ + Ek — Eqct — Eq-o^ — U 



^ — (-k — + E'o-cr 

= 



LO + ek — EQa — Eo-a — U 



Now one can substitute the GF in ( p9lD by the expression ( 299] ). This gives 
a new approximate dynamic solution of SIAM where the complex expansion 
both in U and V was incorporated. The important observation is that this 
new solution satisfies both the limits ( |272| ). For example, if we wish to get a 
lowest order approximation up to C/^ and V"^, it is very easy to notice that 
for V = 0: 



« foacl^^Ck-alfoa^O-a 
(304) « Cfc_^4_^/oa|/olno-a 



< c\_^Ck-a X no-a > 



UJ 



E, 



Oct 



U 



ijj 



Eoa-U 



This results in the possibility to find explicitly all necessary quantities and, 
thus, to solve the problem in a self-consistent way. 

In summary, we presented here a consistent many-body approach to analytic 
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dynamic solution of SIAM at finite temperatures and for a broad interval 
of the values of the model parameters. We used the exact result ( |289| ) to 
connect the single-particle GF with higher-order GF to obtain a complex 
combined expansion in terms of U and V for the propagator. To summa- 
rize, we reformulated the problem of searches for an appropriate many-body 
dynamic solution for SIAM in a way that provides us with an effective and 
workable scheme for constructing of advanced analytic approximate solu- 
tions for the single-particle GFs on the level of the higher-order GFs in a 
rather systematic self-consistent way. This procedure has the advantage 
that it systematically uses the principle of interpolation solution within the 
equation-of-motion approach for GFs. The leading principle, which we used 
here was to look more carefully for the intrinsic functional structure of the 
required relevant solution and then to formulate approximations for the 
higher-order GFs in accordance with this structure. 

The main results of our IGF study are the exact Dyson equations for the 
full matrix GFs and a new derivation of the GMF GFs . The approximate 
explicit calculations of inelastic self-energy corrections are quite straightfor- 
ward but tedious and too extended for their description. Here we want to 
emphasize an essentially new point of view on the derivation of the General- 
ized Mean Fields for SIAM when we are interested in the interpolation finite 
temperature solution for the single-particle propagator. Our final solutions 
have the correct functional structure and differ essentially from previous so- 
lutions. 

Of course, there are important criteria to be met (mainly numerically) , such 
as the question left open, whether the present approximation satisfies the 
Friedel sum rule (this question left open in |105| ] and [ |107[ | ) . A quantitative 
numerical comparison of self-consistent results {e.g. the width and shape of 
the Kondo resonance in the near-integer regime of the SIAM) would be cru- 
cial too. In the present consideration, we concentrated on the problem of 
correct functional structure of the single-particle GF itself. In addition to 
SIAM, it will be instructive to consider sketchy the PAM and TIAM too for 
completness. 

9.8 Quasi-Particle Dynamics of PAM 

The main drawback of the H-F type solution of PAM ( pl| ) is that it ignores 
the correlations of the ""up" and "down" electrons. In this Section, we take 
into account the latter correlations in a self-consistent way using the IGF 
method. We consider the relevant matrix GF of the form ( cf. ( p24D ) 



(305) 
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The equation of motion for GF (305) reads 



-Vk 



«Cfc,|cL» «Cka\fL» 



(306) 



-Vk (W -Ek)J\ « /fc^|4^ » « fka\fL » ) 



«.4|4» «A\fl» 



^ pq 



where A = fk+pafl+q-afq-t^- According to eq.(30), the definition of the 
irreducible parts in the equation of motion ( 306D are defined as fohows 

^ fk+pafl+q-afq-crl'^kij >> = << fk+pafp+q-afq-A^ka ~ 

5p,Q < nq^„ >« fka\c\^» 

fk+pafp+q-afq-Afk(j >> = << fk+pafp+q-afq-Afka ~ 

< ^q-o >« fkalfka » 

After substituting these definitions into equation (p06D, we obtain 



- efc) 



-Vk 



' « cfcalcL » « ckAfL » 

Vk {^-EAk))J \« fka\cl» «fka\fL», 



(307) 



pq 





A\r\ >> (*^) ^i^"!" 



« A\cl^ » « A\fl » 



The following notation was used 

E„{k) = Ek- Un_„; U-a =< fl-afk-a > 

The definition of the generalized mean field GF ( which, for the weak 
Coulomb correlation [/, coincides with the Hartree-Fock mean field ) is evi- 
dent. All inelastic renormalization terms are now related to the last term in 
the equation of motion (|307] ). All elastic scattering ( or mean field) renor- 
malization terms are included into the following mean-field GF 



[y^ - ek) 
-Vk 



-Vk \(« cfc<xicL »° « ckAfL »° 



1 
1 



(308) 

It is easy to find that (c/. ( |229D and (p30|)) 
(309) « fkAfl »°= - E^{k) 



\Vk\' A-i 

UJ-Ck 

\Vk\' 
uJ-E,{k) 

At this point, it is worthwhile to emphasize a significant difference between 
both the models, PAM and SIAM. The corresponding SIAM equation for 



(310) « Cfc<,|4^ »°= Iw-efc 
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generalized mean field GF (228) reads 



E 



•(a; 



(311) 



-Vp5pk ^ « ^ka\cl^ »° « Cka\fL 

Eo^ - Un.^) ][« /o,|cL »° « /o.l4. »° 



This matrix notation for SIAM shows a fundamental distinction between 
SIAM and PAM. For SIAM, we have a different number of states for a 
strongly localized level and the conduction electron subsystem: the conduc- 
tion band contains 2A^ states, whereas the localized (s-type) level contains 

shows clearly that this diffi- 



only two. The comparison of ( 311 ) and ( |30^ 
culty does not exist for PAM : the number of states both in the localized 
and itinerant subsystems are the same, i.e. 2N. 

This important difference between SIAM and PAM appears also when we 
calculate inelastic scattering or self-energy corrections. By analogy with the 
Hubbard model, the equation of motion (307) for PAM can be transformed 
exactly to the scattering equation of the form ( 36). Then, we are able to 
write down explicitly the Dyson equation (37) and the exact expression for 
the self-energy M in the matrix form: 





M22 



(312) MM -- 

Here the matrix element M22 is of the form 



(313) 



22 



El 

iV2 



''<< fk-\-pafl+q-afq-o-\fl-afr+s-af) 



o" J k+sa 



pqrs 



To calculate the self-energy operator ( pi 3D in a self-consistent way, we pro- 
ceed by analogy with the Hubbard model in Section 8.1. Then we find both 
expressions for the self-energy operator in form ( |149| ) and ( [1521) . 



1 0^ 
1 



9.9 Quasi-Particle Dynamics of TIAM 

Let us see now how to rewrite the results of the preceeding Sections for the 
case of TIAM Hamiltonian (62). We again consider the relevant matrix GF 
of the form (cf.(|224D ) 

/ « Cfc^lcL » « Cfcal/L » « Cka\fL » 

= « haWka » « h^lfL » « flalfL » 
\ « /2a|cL » « f2a\fL » « f2a\fL » , 
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The equation of motion for GF (314) reads 



E 

p 

(315) 



■ (w - ep)6pk 



-Vip 

-V2p 



V2I 



— Vlp^pk 
-V12 




+ 



«^2|cL» «^2|/L» «^2|/L»- 

The notation is as follows 

M = flafl-afl-cr'i M = f2af2-ah~a 

In a compact notation, the equation (|315D has the form (c/. ( p55|) ) 
(316) ^F(p,fe)Gpfc(a;) = / + C/Dp(a;) 



We thus have the equatin of motion ( |316| ) which is a complete analogue of 
the corresponding equations for the SIAM and PAM. After introducing the 
irreducible parts by analogy with the equation (225) 

- < ><< fi„\B »^ 

- < n2-a >« f2a\B >><^ 



and performing the second-time differentiation of the higher-order GF, and 
introducing the relevant irreducible parts, the equation of motion ( |316| ) is 
rewritten in the form of Dyson equation (37). The definition of the gener- 
alized mean field GF is as follows 



E 

p 

(317) 



-Vip6pk 



-Vip 

-V2p 



-V21 



-VipSpk 
-V12 










<x) 


<-^12 


Gl3\ I 


n 





<-^22 




1 





<-^32 


^^3/ ^ 


VO 


1 



The matrix GF (317) describes the mean-field solution of the TIAM Hamil- 
tonian. The explicit solutions for diagonal elements of are ( cf. ( p29[ )) 



(318);< Cka\cl^ »^= (t^ - efc 



\Vik\' 



to - [Eq^ - Un-a) 



Au{k,uj) 



-1 
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(319) « fiMl »°= - (Eoa - Un^„) - S{u)) - A22(fe,u;)) 

(320) « /2.I/2I »°= (uJ - {Eoa ~ Un^^) ~ S{lo)) - AgalA:, u;)) 
Here we introduced the notation 



V21V12 ,_1 



\y |2 

A22(A:,u;) = {\2i{oo) + Fi2)(A2i(u;) + V2i)[uo - {E^, - Un.,) - 

io — ep 

A33{k,uj) = (Ai2(a;) + F2i)(Ai2(a;) + Vi2)[uJ - {Eoa - Un^^) 



Epl^ipi 



(321) A12 = A21 = ^ 



The formal solution of the Dyson equation for TIAM contains the self-energy 
matrix 

/O 

(322) M = p M22 M23 

\0 M32 M33, 

where 

(323) M22 = [/'((''-)« /i.ni_.|//^ni_. »(-))^' 
M32 = [/^((^'■)« /2<xn2-.|/Lni_, 

M23 = \J\^'''-K< han^-a\fln2-. »^''^Y 

M33 = [/^((^^•)« /2<xn2_<x|/2U2-a 

To calculate the matrix elements ( p23D, t he same procedure can be used as 
it was done previously for the SIAM ( p39|) . As a result, we find the following 
explicit expressions for the self-energy matrix elements ( cf.(241) 

J-00 UO — UJi - UJ2 

{--Im « S{\Sl >><^i) 
vr 



(324) « /^J/t^ 

M22(w) = U CLjJidAjJ2 
J —00 



-00 UJ — Wi — LJ2 

{--Im « Si\S{ »u 
vr 
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1 



(325) (--/m«/it|/i\ » 



TT 



J— oo 



1 + iV(a;i) - n(u;2) 



oo a; — wi — a;2 



(326) (__/^ « /,J/T^ »^^^ 



TT 



2 i+^('^i)-"('^2) 



/■+00 

Afayw) = C/^ / (ia;i(ia;2 

J— oo 



a; — — 6<J2 

TT 

(327) |/t»^j 

TT ' 

where the following notation was used: 



For M33 we obtain the same expressions as for M22 with the substitution 
of index 1 by 2. For Mgj we must do the same. It is possible to say that 
the diagonal elements M22 and M33 describe single-site inelastic scattering 
processes; off-diagonal elements M23 and M32 describe intersite inelastic 
scattering processes. They are responsible for the specific features of the 
dynamic behaviour of TIAM ( as well as the off-diagonal matrix elements 
of the GF G^) and, more generally, the cluster impurity Anderson model 
(CIAM). The nonlocal contributions to the total spin susceptibility of two 
well formed impurity magnetic moments at a distance R can be estimated 
as 

(328) xpair -« ST\St »~ 2x - 127rS^(^)^^^^ 

giJ,B [kFRr 

In the region of interplay of the RKKY and Kondo behaviour, the key point 
is then to connect the partial Kondo screening effects with the low tempera- 
ture behaviour of the total spin susceptibility. As it is known, it is quite dif- 
ficult to describe such a threshold behaviour analytically. However, progress 
is expected due to a better understanding of the quasi-particle many-body 
dynamics both from analytical and numerical investigations. 



10 Conclusions 

In the present paper, we have formulated the theory of the correlation effects 
for many-particle interacting systems using the ideas of quantum field the- 
ory for interacting electron and spin systems on a lattice. The workable and 
self-consistent IGF approach to the decoupling problem for the cquation- 
of-motion method for double-time temperature Green functions has been 
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presented. The main achievement of this formulation was the derivation 
of the Dyson equation for double-time retarded Green functions instead of 
causal ones. That formulation permits to unify convenient analytical prop- 
erties of retarded and advanced GF and the formal solution of the Dyson 
equation (38), that, in spite of the required approximations for the self- 
energy, provides the correct functional structure of single-particle GF. The 
main advantage of the mathematical formalism is brought out by show- 
ing how elastic scattering corrections (generalized mean fields) and inelastic 
scattering effects (damping and finite lifetimes) could be self-consistently 
incorporated in a general and compact manner. In this paper, we have 
thoroughly considered the idealized Anderson and Hubbard models which 
are the simplest (in the sense of formulation, but not solution) and most 
popular models of correlated lattice fermions. We have presented here the 
novel method of calculation of quasi-particle spectra for these and basic spin 
lattice models, as the most representative examples. Using the IGF method, 
we were able to obtain a closed self-consistent set of equations determining 
the electron GF and self-energy. For the Hubbard and Anderson models, 
these equations give a general microscopic description of correlation effects 
both for the weak and strong Coulomb correlation, and, thus, determine the 
interpolation solutions of the models. Moreover, this approach gives the 
workable scheme for the definition of relevant generalized mean fields writ- 
ten in terms of appropriate correlators. 

We hope that these considerations have been done with sufficient details 
to bring out their scope and power, since we believe that this technique 
will have application to a variety of many-body systems with complicated 
spectra and strong interaction. The application of the IGF method to the 
investigation of nonlocal correlations and quasi-particle interactions in An- 
derson models has a particular interest for studying of the intersite 
correlation effects in the concentrated Kondo system . A comparative study 
of real many-body dynamics of single-impurity, two-impurity, and periodic 
Anderson model, especially for strong but finite Coulomb correlation, when 
perturbation expansion in U does not work, is important to characterize the 
true quasi-particle excitations and the role of magnetic correlations. It was 
shown that the physics of two-impurity Anderson model can be understood 
in terms of competition between itinerant motion of carriers and magnetic 
correlations of the RKKY nature. This issue is still very controversial and 
the additional efforts must be applied in this field. 

The application of the IGF method to the theory of magnetic semiconduc- 
tors was done in [^], |27|. As a remarkable result of our approach, let us 
mention the generalization of the Shastry-Mattis theory for the magnetic 
polaron to the finite temperatures [p^ . The quasi-particle many-body dy- 



namics of ferromagnetic 1 26] and aniferromagnetic semiconductors [ |l09| , | IIC ] 
was studied too. These studies clarified greatly the true nature of carriers 
in magnetic semiconductors. The application of the IGF method to gener- 
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alized spin-fermion models that was made in papers [^],[111| allows one to 
consider carefully the true nature of carriers in oxides and rare-earth metals. 
These applications illustrate some of subtle details of the IGF approach and 
exhibit their physical significance in a representative form. 
As it is seen, this treatment has advantages in comparison with the standard 
methods of decoupling of higher-order GFs within the equation-of-motion 
approach, namely, the following: 

(i) At the mean-field level, the GF, one obtains, is richer than that fol- 
lowing from the standard procedures. The generalized mean fields 
represent all elastic scattering renormalizations in a compact form. 

(ii) The approximations ( the decoupling ) are introduced at a later stage 
with respect to other methods, i.e. only into the rigorously obtained 
self-energy. 

(iii) Many standard results of the many-particle system theory are repro- 
duced mathematically incomparable more simply. 

(iv) The physical picture of elastic and inelastic scattering processes in 
the interacting many-particle systems is clearly seen at every stage 
of calculations, which is not the case with the standard methods of 
decoupling. 

(v) The main advantage of the whole method is the possibility of a self- 
consistent description of quasi-particle spectra and their damping in a 
unified and coherent fashion. 

(vi) This new picture of interacting many-particle systems on a lattice is 
far richer and gives more possibilities for the analysis of phenomena 
which can actually take place. In this sense the approach we suggest 
produces more advanced physical picture of the quasi-particle many- 
body dynamics. 



Despite the novelty of the IGF techniques introduced above and some (not 
really big) complexity of the details in its demonstrations, the major conclu- 
sions of the present paper can be made intelligible to any reader. The most 
important conclusion to be drawn from the present consideration is that 
the GMF for the case of strong Coulomb interaction has quite a nontrivial 
structure and cannot be reduced to the mean-density functional. This last 
statement resembles very much the situation with strongly non-equilibrium 
systems, where only the single-particle distribution function is insufficient 
to describe the essence of the strongly non-equilibrium state Therefore a 
more complicated correlation functions are to be taken into account, in ac- 
cordance with general ideas of Bogoliubov and Mori-Zwanzig. The IGF 



92 



method is intimately related to the projection method in the sense, that it 
expresses the idea of "reduced description" of a system in the most general 
form. This line of consideration is very promising for developing the com- 
plete and self-contained theory of strongly interacting many-body systems 
on a lattice. Our main results reveal the fundamental importance of the 
adequate definition of Generalized Mean Fields at finite temperatures, that 
results in a deeper insight into the nature of quasi-particle states of the cor- 
related lattice fermions and spins. We believe that our approach offers a 
new way for systematic constructions of the approximate dynamic solutions 
of the Hubbard, SIAM, TIAM, PAM, spin-fermion, and other models of the 
strongly correlated electron systems on a lattice. The work in this direction 
is in progress. 
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A Appendix . The Gram-Schmidt Orthogonalization Pro- 



In this appendix we briefly recall the Gram-Schmidt Orthogonalization Pro- 
cedure. The Gram-Schmidt orthogonalization procedure is an inductive 
technique to generate a mutually orthogonal set from any linearly inde- 
pendent set of vectors. 

Suppose we have an arbitrary n-dimensional Euclidean space, which means 
that scalar multiplication has been introduced in some fashion into an n- 
dimensional linear space. The vectors / and g are orthogonal if their scalar 
product is zero 



We now describe the orthogonalization process, which is a means of passing 
from any linearly independent system of k vectors /i, /2, ...fk to an orthog- 
onal system, also consisting of k nonzero vectors. We denote these vectors 
by 31,52, -5fe- 

Let us put gi = fi, which is to say that the first vector of our system will 
enter into the orthogonal system we are building. After that, put 



cedure 



(Ai) 



{A.2) 



g2 = h + "51 
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Since gi = fi and the vectors /i and /2 are linearly independent, it follows 
that the vector §2 is different from zero for any scalar a. We choose this 
scalar from the constraint 



(A3) = (51,52) = a(5i,5i) + (5i,/2) 

whence 

(91,12) 



(A.4) a = - 



{91,91] 



In other words, we get 52 by subtracting from /2 the projection of /2 onto 
gi. Proceeding inductively, we find 

(A5) = 

j^i \9j,9j) 

We are left with mutually orthogonal vectors which have the same span as 
the original set. 

Let us consider an important example of a basis /i, /2, /s, /4 in a 4-dimensional 
space and then construct the orthonormal basis of the same space. Next, in 
the equality 53 = /s + Pigi + f3292, chose (3i and P2 such that the conditions 
93^91,93^92 are fulfilled. 
From the equalities 

(A6) (51,53) = {91, f3)+Pi{9i, 91) + 02(91,92) 



(A.7) (92,93) = (92, f 3) + 1^1(91,92) + 1^2(92,92) 

we obtain 

(A8) p, = J9ij3l . (92, h) 



(91,91) (92,92) 
Finally, from the equality gi = fi+ 7151 + 7252 + 7353 we find 

/.ON (91, fi) (92, fi) (93, f 4) 
[A.9) 71 = ; 72 = r ; 73 = -7 r 

(51,51) (52,52) (53,53) 

Thus, we see that with the choice of q;,/3i,/32,7i,72,73 made, the vectors 
51,52,53,54 are pairwise orthogonal. 



94 



B Appendix . Moments and Green Functions 



It is known that the method of moments [112] of spectral density is con- 
sidered sometimes as an alternative approach for describing the many-body 
quasi-particle dynamics of interacting many-particle systems. The moments 
technique appears naturally when studying the particle dynamics in many- 
particle systems in the context of time-dependent correlation functions ( 
magnetic resonance, liquids, etc.). Qualitatively, a correlation function de- 
scribes how long a given property of a system persists until it is averaged out 
by the microscopic motion of particles in the macroscopic system. The time 
dependence of a particle correlation function sometimes is approximated (at 
small times) via a power series expansion about the initial time 0. 

cxD _ifi Jn 

{B.l) < Amit) >= ^ -— < A{0)A{t) > \t=o = 

n=0 

Y.^-^ < AmHAH...[H.Am..]]\ > 

Oil. 

The spectral theorem (26), (27) connects A{uj) and the correlation functions. 
From the above expression we obtain the moments M„ of the spectral density 
function 

{B.2) Mn = ^ r dcvcv'^Aiu;) = (-1)" < [[H, [H...[H, A]...]], B]r, > 

So, by definition, the moments are time-independent correlation functions of 
a combination of the operators. In principle, it is possible to calculate them 
in a regular way; however, in practice, it is possible to do this only for a first 
few moments. If the moments of a given spectral density form a positive 
sequence, then GF of appropriate operators is a limit of the sequence 

(i?.3) G{E)= Vim^Gn{E,^) 

Here the parameter — oo < 7 < -|-oo and is real. The approximation pro- 
cedure for GF consists in replacing the G{E) by Gn{E,j), that depends 
also on the appropriate choice of the parameter 7. The Gn{E,^) have the 
properties 

{BA) Gn{E,^) = Gn^i{E,0) 

and are represented by the fraction 

Qn+l{E) - jQn{E) 



(B.b) Gn{E,^)=Mo 



Pn+l{E) - -iPn{E) 
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The polynomials P„ are given by the determinant 



(B.6) 



Pn>liE) 



Mq Ml 
Ml M2 



Mn-l 
1 



Mn 

E 



Mn+1 

M-2n-l 



where 



{B.7) 



D, 



n>l 



Po = l 



Mo Ml 
Ml M2 



Mn 

Mn+l 



Mn Mn+l ... M2n 

L»o = -D-1 = Mo 

The polynomial Qn{E) (which is of (n-l)-th order in E) is related to the 
polynomial Pn{E) (which is of n-th order in E) via the following relation 

Pn{E)-Pn{uj) 



iB.8) 



Qn{E) 



1 r 

tM^J-, 



27rMo J -00 E-LO 
It is possible to find a few lowest-order terms 



-A{Lo)du 



PoiE) = 1; Pi{E) = 



E-Mi 

^ Mq 
M2 - Mq-^ 



(5.10) 



Qo{E) = 0; Qi = 



M2 - Mq-i 

The expression (B.5) can be represented in the following form 



(B.ll) 



G„(^,,)=MoE ^^^^^ 



t{E-Ei{j) 
Here the numbers Ei{^) are roots of the equation 

(B.12) Pn+l{E)-^Pn{E)=Q 

These relations lead to the possibility of practical applications of the moment 
expansion method. If we know the first (2n + 2) moments, then the equation 
(B.12) determines (n + 1) different roots Ei('y). Thus, the spectral density 
function can be represented by 



n+l 



(S.13) 



A{uj) = 27rMo ^ miS{uj - Ei) 



i=l 
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For example, if we know the moments Mq, Mi, M2 then we find, from the 
equation (B.ll), the roots of (B.12) 

{B.U) ^1(7) = MiMq-i + 7(M2 - Mg-i) 

In this approximation, the GF and corresponding spectral density are rep- 
resented as 

(5.15) Go{7) = ^^^^; A{u:) = 2^M^6{iv-Ei) 

It is clear that the Tyablikov decoupling approximation (43) corresponds to 
this approximation within the moment method. An improved decoupling 
scheme, that conserves the first several frequency moments of the spectral 
weight function for the Heisenberg and Hubbard models was developed in 
paper|l| (c/. 

It was shown in ref. ||2q| that the IGF method permits one to calculate 
the spectral density for the spin-fermion model in the approximation that 
preserves the first four moments. This is valid also for the approximation 
used for the strongly correlated Hubbard model in Section 7.2. 
It must be clear from the above consideration that the structure of the 
obtained solution for single-particle GF depends strongly on the stage at 
which irreducible parts were introduced |25|. To clarify this, let us consider 
equation (29) again. Instead of (30), we introduce now the IGFs in the 
following way 

ujG{oj) = Mo+ « [A,H]^ I »^ 

(BAG) uj « [A,H]\A^ »= Ml + (("')« [[A,H]H] \ A^ »^) + 

ai « A\A^ »^ +a2 « [A,H]\A'' »^ 

The unknown constants ai and 02 are connected by the orthogonality con- 
dition 

{BA7) <[[[A,H]H]^''~\A^]>=0 

For illustration, we consider the simplest possibility and write down the 
following equation 

{BA8) w((^'")« [[A,H]H]\A^ ») = ((^'■)« [[A, H]H]\[H, A^] ») 

Then by introducing the irreducible parts for the right operators, we obtain 
{BA9) 

((^")« [[A,H]H]\A^ »){io-a\) = ((^")« [[A, H]H]\[H, A^] »(^")) 
It is clear enough that, as a result, we arrive at the following set of equations 
uj « AIA'' »^ - « [A,H]_ I A'f »uj= Mo 



97 



{B.20) ai «A\A^ »^ +{uj - 02) « [A,H]\A^ Mi - $ 

where 

(5.21) $ = ((*^)« [[A,H]H] I [A,H]^ 

The solutions of the equations (B.20) are given by 
(B.22) « AlA* »„= M„l.-a.)- {Ah-<l» 

{iJ.23) « [^.gl|At»„= "'*f'-'^> + ''-"° 

UJ{UJ — a2) + «! 

aiMo + Q2M1 = M2 

It is evident that there is similarity between the obtained solutions and the 
moment expansion method. The structure of equation (B.22) corresponds 
to the moment expansion (B.ll) except for the factor <I> that should be 
calculated by considering high-order equations of motion or by some relevant 
approximation. 

C Appendix . Projection methods and IGFs 

The IGFs method is intimately related to the projection operator method 
[50 1, [^], that incorporates the idea of "reduced description" of a system 



in the most suitable form. The projection operation [114], [|51[| makes it 
possible to reduce the infinite hierarchy of coupled equations to a few rel- 
atively simple equations that "effectively" take into account the essential 
information about the system that determines the specific nature of the 
given problem. Projection techniques become standard in the study of cer- 
tain dynamic processes. Projection operator techniques of Mori-Zwanzig 



and similar ones [48| are useful for the derivation of relaxation equations 
and formulas for transport coefficients in terms of microscopic properties. 
This approach was applied to a large variety of phenomena concerning the 
line-shape problem. It was shown that there is a close relationship between 
the Mori procedure and the "classical moment problem" of mathematical 
analysis. 

Let us briefly consider the projection formalism for double-time retarded 



GFs [114|, [51|. Ichiyanagi [114| constructed the following set of equations 

for OF (28): 

(C.l) 

(^ - iuk) « Akit),Al{t') »= -i5{t - t') < [Ak,Al] > +F{k, t - t') 
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{C.2) (- + iuok)F{k, t-t') = +i6{t - t') < [K{k), 4] > +n(A;, t - t') 
at 

where F{k, t-t') =« K{k, t),Al{t') » and U{k, t-t') =« K{k, t), K^{k, t') » 
Here, the definitions were introduced: 

(C.3) i^, = SSMhl^^^ Kik,t) = {l-P)A,it) 



{CA) PG =< [G, 4] X [Ak, 4] Ak 

The projection operator P defined in (C.4) is different from the one intro- 
duced by Mori. The main result of paper []114|| is that, using the projection 
operator, a Dyson equation that determines an irreducible quantity, proper 
self-energy part, was obtained in the following form: 

(as) {u;-u;k- —^^^-r—M{k,u)) « A,\aI »^= _ < i^t^A] > 

<[^fc,4]> 2^ 

Here M{k,oj) is the self-energy , that, in the diagrammatic language, con- 
sists of irreducible diagrams. 



Our point of view is closely related to that of ref. [114| and to the develop- 
ment of the this paper by Tserkovnikov in a systematic way |51]. However, 
our strategy is slightly different in the time evolution aspect. We consider 
our IGF technique as more convenient from the practical computational 
point of view. 

D Appendix . Effective Perturbation Expansion for the 
Mass Operator 

Let us consider a useful example how to iterate the initial "trial" solution 



and to get an expansion for the mass operator To be concrete, let 
us consider the calculation of the mass operator for the Hubbard model in 
Section 8.1. The first iteration for the equation ( |147D with the trial function 



(|148| ) have lead us to the expression (149), which we rewrite here in the 
following form 



pq 

where 



■^kpq — ''^p+g— (t(1 n-k+pa 1^q—a) ~l" ^k+pa 



n. 



q-a 
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^kpq = -e(p + qa) + e{k + pa) + e{qa) 
Now we are able to calculate the spectral weight function gka{'^) (0) 

We approximate this expression by the following way 
(Z?.3) gUu^) ^{l- ak.)5{u; - E^,) + -- ^"""^"^^ 



Here 

r, r,,,"! — ^ 

Eku = e{ka) + Aka 
Aka = ReMkai'^ + ie) 
The unknown factor (1 — aka) is determined by the normalization condition 



f/2 

rfccr(w) = ^ Nkpq6{uJ - ^kpq) 

pq 



/oo 
-oo 



whence 



_ ^ \ ^ Hpq 



Then, using (22), we find for the mean occupation numbers 

Now we can use the spectral weight function ( D.2) to iterate the equation 
( |147] ) and to get a perturbation expansion for the self-energy M^a in the 
pair approximation. Instead of the initial trial solution in the form of delta- 
function ( |148[) , we take the expression (D.3). It is easy to check that we get 
an expansion up to 6th order in U. 
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Table 1: EvolutioTi of the Mean Field Concept. 



Type of the mean field 

Uniform molecular field 

in dense gases 

Uniform internal mean field 

in magnets 

Thomas-Fermi model 

Uniform mean field 

in many-electron atoms 

Molecular mean field in 

Heisenberg ferromagnet 

Non-uniform(local) staggered 

mean field in antiferromagnet 

Reaction and cavity field in 

polar substances 

Stoner mean-field model 

of band magnetism 

Slater mean-field model 

of band antiferromagnet ism 

BCS-Bogoliubov mean field in 

superconductors 

Tyablikov decoupling 

for Heisenberg ferromagnet 

Mean field theory for SIAM 

Density Functional Theory 

for inhomogeneous electron gas 

C alien decoupling 

for Heisenberg ferromagnet 

Alloy analogy (mean field) approximation 

in strongly correlated model 

Generalized Mean Fields 

in Heisenberg ferromagnet 

Spin Glass Mean Field Model 

Generalised Mean Fields in Strongly 

Correlated Hubbard Model 

Generalized Mean Fields 

in Heisenberg antiferromagnet 

Generalized Mean Fields 

for itinerant antiferromagnet 



Author Year 

van der Waals 1873 

P.Weiss 1905 

L.H.Thomas, E.Fermi 1926-28 

D.Hartree , V.Fock 1928-32 

F.Bloch 1930 

L.Neel 1932 

L.Onsager 1936 

E.Stoner 1938 

J.Slater 1951 

N.N.Bogoliubov 1958 

S. Tyablikov 1959 

P.W.Anderson 1961 

W.Kohn 1964 

H.B.Callen 1964 

J.Hubbard 1964 

N.Plakida 1973 

S.F.Edwards, P.W.Anderson 1975 

A.L.Kuzemsky 1975-78 

A.L.Kuzemsky, D.Marvakov 1990 

A.L.Kuzemsky 1999 
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